Draft version January 20, 2012 

Preprint typeset using L2TpX style emulateapj v. 5/2/1 1 



THE HUBBLE SPACE TELESCOPE" CLUSTER SUPERNOVA SURVEY: 
III. CORRELATED PROPERTIES OF TYPE IA SUPERNOVAE AND THEIR HOSTS AT 0.9 < z < 1.46 

J. Meyers 1,2 , G. Aldering 2 , K. Barbary 1 - 2 , L. F. Barrientos 3 , M. Brodwin 4 - 21 , K. S. Dawson 5 , S. Deustua 6 , M. Doi 7 , 
P. Eisenhardt 8 , L. Faccioli 2 , H. K. Fakhouri 1 ' 2 , A. S. Fruchter 6 , D. G. Gilbank 9 , M. D. Gladders 10 , G. GOLDHABER 1 ' 2 ' 23 , 
A. H. Gonzalez", T. Hattori 12 , E. Hsiao 2 , Y. Ihara 7 - 22 , N. Kashikawa 13 , B. Koester 1014 , K. Konishi 15 , C. Lidman 16 , 

L. LUBIN 17 , T. MOROKUMA 713 ' 22 , T. ODA 18 , S. PERLMUTTEr'' 2 , M. POSTMAN 6 , P. RlPOCHE 2 , P. ROSATI 19 , D. RUBIN 1 ' 2 , E. RYKOFF 2 , 

A. Spadafora 2 , S. A. Stanford 17 20 , N. Suzuki 2 , N. Takanashi 13 , K. Tokita 7 , N. Yasuda 15 
(The Supernova Cosmology Project) 

Draft version January 20, 2012 

ABSTRACT 

Using the sample of Type la supernovae (SNe la) discovered by the Hubble Space Telescope (HST) Cluster 
Supernova Survey and augmented with ffi>r-observed SNe la in the GOODS fields, we search for correlations 
between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative 
morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving 
early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify 
six SN la hosts that are early-type cluster members and eleven SN la hosts that are early-type field galaxies. 
We confirm for the first time at z > 0.9 that SNe la hosted by early-type galaxies brighten and fade more 
quickly than SNe la hosted by late-type galaxies. We also show that the two samples of hosts produce SNe la 
with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive 
galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates 
of late-type GOODS SN hosts from Thomson & Chary (201 1), we investigate the correlation of host mass with 
Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint 
of this relation is found at z > 0.9. By simultaneously fitting the average cluster galaxy formation history and 
dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely 
E(B — V) < 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type 
galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted 
SNe la studied here occupy a more favorable environment to use as well-characterized high-redshift standard 
candles than other SNe la. 

Subject headings: Cosmology:General, Supernovae:General, Galaxies:Clusters 
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1. INTRODUCTION 

The use of Type la supernovae (SNe la) as standard candles 
in estimating astronomical distances has proven indispensable 
for modern cosmology, leading to the remarkable discovery 
that the expansion of the Universe is accelerating (Riess et al. 
1998; Perlmutter et al. 1999). Together with measurements of 
the cosmic microwave background (Dunkley et al. 2009; Ko- 
matsu et al. 201 1) and baryon acoustic oscillations (Eisenstein 
et al. 2005; Percival et al. 2010), our recent compilation of the 
myriad SN surveys that led to and followed the initial discov- 
ery has now constrained the dark energy equation of state to 
bew= -0.974t°;p^(stat)t°;°^(stat+sys) (Amanullah et al. 
2010). The constraints from the 557 SNe la in this compila- 
tion are now limited by systematic uncertainties. Additional 
SN la observations will not greatly improve dark energy con- 
straints until the systematic errors are reduced. 

Improvements in SN la distance estimates have historically 
followed two paths. The first path is to identify subsets of SNe 
with smaller intrinsic brightness dispersion (e.g. the Type la 
subset or "Branch-normal" SNe la). The second path is to 
make corrections exploiting empirical correlations between 
SN brightness and other observables (e.g. light curve shape 
and color). The properties of SN la host galaxies offer ad- 
ditional variables with which to search for smaller dispersion 
subsets or correlated observables. In fact, such a correlation 
has recently been reported: SNe la in more massive hosts are 
brighter after light curve shape and color corrections (Kelly 
et al. 2010; Sullivan et al. 2010; Lampeitl et al. 2010). 

The brightnesses of SNe la are correlated with their colors 
at maximum light: bluer SNe la are brighter. This trend is 
broadly consistent with extinction due to dust as commonly 
described by R B = A{B)/E{B - V) (Cardelli et al. 1989). 
When relating SN la colors to peak magnitudes, the param- 
eter /3 = AMb/Ac is used, where c = (B-V) - 
(B — V) max ; e.g. c is the excess color at maximum bright- 
ness with respect to the average. Direct fits of SN la observa- 
tions have found that j3 is considerably smaller than the value 
of Rb observed in the Milky Way diffuse interstellar medium, 
suggesting that an intrinsic component to the bluer-brighter 
relation is also important (Tripp 1998; Astier et al. 2006; Con- 
ley et al. 2007; Nobili & Goobar 2008). The systematic errors 
associated with the interpretation of j3, i.e., what fraction is 
attributable to dust and how might that evolve with redshift, 
are now comparable to the statistical errors of SN la surveys 
(Kowalski et al. 2008; Wood-Vasey et al. 2007). However, if 
one were to identify SN la hosts with small extinction, then a 
single component of the SN la color-magnitude relation, the 
intrinsic component, could be isolated. 

In this paper we analyze the host galaxies of z > 0.9 SNe la 
drawn from the HST Cluster SN Survey (Dawson et al. 2009, 
Pl-Perlmutter: GO- 10496) and surveys of the GOODS fields 
(Blakeslee et al. 2003b; Riess et al. 2004, 2007). We investi- 
gate how host properties correlate with SN la properties and 
identify hosts minimally affected by dust. Particularly inter- 
esting host properties such as age and metallicity are usually 
difficult to obtain, but are correlated with the position of a 
galaxy on its cluster's red sequence (Gallazzi et al. 2006). 
The environment of the HST Cluster SN Survey thus lends 
itself to a particular analysis strategy: how does the location 
of a host galaxy on its cluster's red sequence correlate with 
its SN, what does it imply about the SN type, and how does it 
begin to parse the origin of the SN reddening law? 

This paper is organized as follows: In §2 we describe 



progress to date of SN host galaxy studies at lower redshifts. 
In §3 we describe the HST Cluster SN Survey and the HST 
SN surveys of the GOODS fields, the data taken, and the var- 
ious photometric, morphological, and spectroscopic quanti- 
ties that we derive from them. In §4 we describe how, by 
using the color-magnitude locations of spectroscopically con- 
firmed red-sequence cluster members, we predict the location 
of the red sequence for each cluster, including those without 
extensive spectroscopic coverage. Additionally, we demon- 
strate how quantitative morphology parameters can be used to 
enhance the contrast of the red sequence to more easily fit a 
color-magnitude relation (CMR) for each cluster. In §5 we de- 
scribe criteria to classify the observed SN hosts, and in §6 we 
use expected rates of SNe la and core collapse SNe (SNe CC) 
in their individual galaxies to help type the SNe. In §7 we 
investigate limits on extinction for SNe hosted by early-type 
galaxies. In §8 we investigate correlations between SN prop- 
erties and properties of their host galaxies. Finally, in §9 we 
summarize our conclusions and discuss their implications. 

This paper is one of a series of ten papers that report 
supernova results from the HST Cluster Supernova Survey 
(PI: Perlmutter, GO- 10496), a survey to discover and follow 
SNe la of very distant clusters. Paper I (Dawson et al. 2009) 
describes the survey strategy and discoveries. Paper II (Bar- 
bary et al. 2010) reports on the SN la rate in clusters. The 
current work, Paper III, addresses the properties of the galax- 
ies that host SNe la. Paper IV Ripoche et al. (2011) intro- 
duces a new technique to calibrate the "zeropoint" of the NIC- 
MOS camera at low count rates, which is critical for placing 
NICMOS-observed SNe la on the Hubble diagram. Paper V 
(Suzuki et al. 201 1) reports the SN la light curves and cosmol- 
ogy from the HST Cluster SN Survey. Paper VI (Barbary et al. 
2011) reports on the volumetric field SN la rate. Melbourne 
et al. (2007), one of several unnumbered papers in this series, 
present a Keck Adaptive Optics observation of a z = 1.31 
SN la in //-band. Barbary et al. (2009) report the discov- 
ery of the extraordinary luminous supernova, SN SCP06F6. 
Morokuma et al. (2010) present the spectroscopic follow-up 
observations for SN candidates. Hsiao et al. (2011) develop 
techniques to remove problematic artifacts remaining after the 
standard STScI pipeline. A separate series of papers, ten to 
date, reports on cluster studies from the HST Cluster SN Sur- 
vey: Brodwin et al. (201 1); Eisenhardt et al. (2008); Jee et al. 
(2009); Hilton et al. (2007, 2009); Huang et al. (2009); Santos 
et al. (2009); Strazzullo et al. (2010); Rosati et al. (2009); and 
Jee etal. (2011). 

2. SNIA HOSTS 

In this section we briefly discuss the relationships between 
SNe and their host galaxies with particular emphasis on early- 
type host galaxies of SNe la. We also introduce the HST 
Cluster SN Survey, which targets massive galaxy clusters to 
increase the yield of SN la discoveries, particularly those in 
low-dust early-type hosts. 

2. 1 . SN typing by host 

Early-type galaxies have almost never been observed to 
host Type Ib/c or Type II SNe (hereafter SNe CC), which are 
generally thought to originate from the collapse of the cores 
of massive (and hence young) stars. In an extensive litera- 
ture search of hundreds of SNe CC, Hakobyan et al. (2008) 
uncovered 22 examples reported as being hosted by early- 
type galaxies, but ultimately reclassified 19 of these galax- 
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ies as late-type. The SN associated with one of the remain- 
ing three galaxies has since been reclassified as a cataclysmic 
variable star within the Milky Way (Leonard 2010). The SNe 
hosted by the remaining two early-type hosts, SN2000ds and 
SN2005cz, are both members of the SN2005E-like subset of 
faint "Ca-rich" SNe lb. The origins of SNe in this subset 
are still under debate, though it is clear that they are read- 
ily distinguishable from SNe la by their faint peak magni- 
tudes (Mb ~ — 15 compared to Mb ~ —18 for even a very 
faint SN la) (Perets et al. 2010; Kawabata et al. 2010). In a 
separate analysis of the near ultraviolet and optical colors of 
early-type SN hosts (including two which were reclassified 
by Hakobyan et al.), Suh et al. (201 1) find that the early-type 
hosts of SNe CC have had more recent star formation and 
occupy a different part of the UV-optical color-magnitude di- 
agram than identically selected early-type hosts of SNe la. 

Despite these rare possible exceptions, SNe hosted by 
early-type galaxies are overwhelmingly Type la. Typing high- 
redshift (z > 0.9) SNe spectroscopically can be done from 
space with HST (Riess et al. 2004) but requires a significant 
investment of orbits (4-8 per SN at z ~ 1.2). Not only are 
SNe increasingly faint at higher redshift, but the rest-frame 
features used to identify a Type la become shifted out of the 
optical wavelengths. Ground-based spectroscopy faces the 
additional difficulty that the night sky becomes increasingly 
bright at redder wavelengths. With ideal conditions and by 
targeting the SN near the peak of its light curve, Morokuma 
et al. (2010) have shown it is possible to spectroscopically 
type SNe la as high as z = 1.34 from the ground, but this 
is by no means typical. In fact, many of the lower redshift 
targets in this paper (of which many also had higher signal- 
to-noise ratio spectra) yielded inconclusive types. Classifying 
SNe in ellipticals galaxies as Type la can therefore provide an 
efficient and robust alternative to potentially expensive spec- 
troscopic typing at high redshift. 

2.2. Dust in early-type SN la hosts 

Early-type SN la hosts offer an additional advantage to typ- 
ing: they generally contain only small amounts of dust. Al- 
though HST images of the cores of nearby early-type galaxies 
indicate that ~ 50% exhibit optical absorption due to dust in 
small disks or filaments, these are usually confined to the cen- 
tral few hundred parsecs (van Dokkum & Franx 1995; Tomita 
et al. 2000; Rest et al. 2001 ; Tran et al. 2001 ; Lauer et al. 2005; 
Ferrarese et al. 2006). In contrast, SNe la hosted by early-type 
galaxies are spatially distributed following the optical light of 
their host galaxies (Forster & Schawinski 2008). Observa- 
tions of emission from dust at far-infrared (FIR) wavelengths 
suggest that additional dust is diffusely distributed through- 
out at least some nearby early-type galaxies (Goudfrooij & de 
Jong 1995; Temi et al. 2004, 2007). The inferred dust mass 
relative to the stellar mass of these galaxies varies greatly; at 
fixed optical luminosity, Temi et al. (2007) measure FIR lumi- 
nosities that span two orders of magnitude. At the low dust- 
mass end, at least several early-type galaxies (3 elliptical and 
3 SO) with M B ~ - 18 to -20, are undetected in the Herschel 
Space Observatory survey of the Virgo cluster. Upper limits 
constrain the associate dust mass of these galaxies to be < 10 4 
M q (Clemens et al. 2010). For comparison, the Virgo clus- 
ter spiral galaxies in the same £>-band luminosity range also 
observed by Herschel have dust masses of ~ 10 7 — 10 8 M ) 
(Davies et al. 2011). Early-type galaxy dust masses derived 
from Spitzer tend to fall in the range 10 4 — 10 6 (Kaneda et al. 
2007). To estimate the extinction associated with a certain 



mass of dust, we use (Ay) — MdT/Y,, with £ the area of the 
dust feature, (Ay) the average U-band absorption in that area, 
and r~6x 10~ 6 mag kpc 2 Mq 1 the visual mass absorption 
coefficient (van Dokkum & Franx 1995). For dust uniformly 
distributed throughout a disc with radius 3 kpc, which is rep- 
resentative of the half-light radii of the galaxies considered in 
this paper, the absorption ranges from Ay < 0.002 to ~ 0.2. 
Since dust is not actually distributed in a uniform foreground 
screen, but rather is embedded within the galaxy, these values 
only provide a rough guide to the expected extinction along 
the lines of sight toward early-type hosted SNe la. 

At high redshift, FIR observations of SN la hosts are im- 
practical. However, optical observations of the integrated 
colors of galaxies can also place constraints on the extinc- 
tion expected from diffuse dust. For example, consider the 
case where galaxy colors are the combined result of stellar 
population age, metallicity, and dust. The contributions of 
age and metallicity to the integrated color can be removed 
by using an auxilliary correlated variable. At low redshift, 
the spectral absorption feature Mg 2 is such a variable. The 
equivalent width of Mg 2 is expected to correlate with the 
rest-frame intrinsic (i.e. unextinguished) B — V color of 
elliptical galaxies (Faber et al. 1989). The line strength in- 
creases with increasing galaxy age and metallicity, and con- 
sequently color, but because the wavelength interval of the 
line is relatively narrow it should not correlate significantly 
with dust. Schlegel et al. (1998) measured the scatter about 
the Mg 2 -color relation for nearby elliptical galaxies to be just 
(Tbv = 0.0257 mag. If at least some elliptical galaxies are 
effectively dust-free (Ay < 0.01), which seems likely given 
the range of dust masses discussed above, then this scatter 
measurement implies an absolute limit on dust reddening at 
the E(B - V) ~ 0.03 level. While spectroscopy of the Mg 2 
feature may not be feasible at high redshift, the red sequence 
may provide a readily available alternative calibrator. 

Two mechanisms have been proposed to explain the origin 
of the dust found in elliptical galaxies. Elliptical galaxies may 
accrete dust and gas during (minor) mergers with dusty late- 
type galaxies. This scenario is supported by observations that 
the motions of dust and gas in some nearby early-type galax- 
ies seem unrelated to the motions of their stars (Goudfrooij 
et al. 1994; van Dokkum & Franx 1995; Caon et al. 2000). 
As part of the Galaxy-Zoo project, Kaviraj et al. (2011) find 
that early-type galaxies with prominent dust lanes are likely 
associated with mergers due to the frequency with which they 
present disturbed morphologies compared to a control sample 
without dust lanes. Prominent dust lane early-type galaxies 
are relatively rare, however, constituting just 4% of the early- 
type Galaxy-Zoo sample, and are even rarer in cluster envi- 
ronments.. 

Noting that the dynamical freefall time for a merging galaxy 
(several ~ 10 8 yr) is comparable to the sputtering lifetime of 
dust grains in the hot interstellar medium of early-type galax- 
ies (~ 10 7 — 10 s yr), Temi et al. (2007) argue that dust orig- 
inating from mergers must be resupplied at least every ~ 10 8 
yr. Such a high rate of early-type - late-type mergers is not 
observed, however. As an alternative, they propose an inter- 
nal origin for diffuse dust in early-type galaxies: dust is gen- 
erated in the atmospheres of evolved red giant stars (Knapp 
et al. 1989; Athey et al. 2002) and may accumulate in the 
concentrated disks of dust commonly observed in the centers 
of early-type galaxies(Mathews & Brighenti 2003). Intermit- 
tent AGN activity can then buoyantly transport this dust out 
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to large radii where it will be destroyed by sputtering in the 
hot interstellar medium. 

One possible caveat to the expectation of low extinction 
SNe la in early-type galaxies is small amounts of recent star- 
formation seen in some early-type galaxies, especially lentic- 
ular galaxies, and associated dust (Yi et al. 2005; Donas 
et al. 2007; Schawinski et al. 2007; Kaviraj et al. 2007, 
2008; Temi et al. 2009a,b). Together with the expected in- 
crease in the SN la rate with star-formation rate (see §2.3 
below), it is possible that an SN la in a mostly dust-free 
galaxy could, nevertheless, suffer from extinction. Fortu- 
nately, recent star-formation in early-type galaxies is usually 
confined to their inner regions (Kuntschner et al. 2006; Sarzi 
et al. 2006), whereas, as mentioned above, the spatial distri- 
bution of SNe la follows the optical light distribution. Fur- 
thermore, we have directly constrained star-formation using 
spectroscopy of our host galaxies and its possible effect in in- 
creasing the galactic SN la rate (see §6). We expect that most 
of our early-type hosted SNe la are not associated with recent 
star-formation. 

2.3. SN la properties by host 

Several lines of evidence suggest that the demographics of 
SNe la hosted by early-type galaxies are different than those 
of SNe la hosted by late-type galaxies. Using infrared magni- 
tudes of low-redshift SN la host galaxies in the 2MASS sur- 
vey, Mannucci et al. (2005) measured the SN la rate per unit 
stellar mass for different host types. They found a rate ~ 20 
times larger in late-type galaxies than in E/S0 galaxies. Split- 
ting the galaxies by color (as a proxy for star formation) they 
found a rate ~ 30 times larger in blue galaxies than in red 
galaxies. Sullivan et al. (2006) fit SED models to five-band 
galaxy photometry to directly constrain the star formation rate 
and mass of intermediate redshift SN la host galaxies in the 
Supernova Legacy Survey (SNLS). This analysis showed that 
the SN la rate is ~ 10 times larger per unit stellar mass in 
actively star-forming galaxies than in passive galaxies. For a 
given episode of star formation, the SN la rate quickly peaks 
and declines as the progenitor population ages. 

Perhaps not surprisingly, the properties of SNe themselves 
are also correlated with the properties of their hosts. At low 
redshift, SNe la hosted by late-type, spiral galaxies tend to 
have broader (slower) light curves on average (Hamuy et al. 
1996; Gallagher et al. 2005), although interestingly, Hicken 
et al. (2009) find that the broadest SN light curves are found in 
SNe hosted by Sb-Sc galaxies and not younger Sd-Irr galax- 
ies. Sullivan et al. (2006) showed that this distinction ex- 
tends to intermediate redshift hosts with actively star-forming 
galaxies (generally late-type) hosting SNe la with broader 
light curves. Several analyses suggest that SNe la hosted by 
early-type galaxies may have a smaller intrinsic peak bright- 
ness dispersion after light curve and color corrections are ap- 
plied (Sullivan et al. 2003; Jha et al. 2007; Sullivan et al. 2010; 
Lampeitl et al. 2010), though at least one analysis finds the 
opposite result (Hicken et al. 2009). 

Since SN la properties are correlated with the properties of 
their host galaxies, and the demographics of galaxies change 
with redshift, the demographics of SNe la will also change 
with redshift. This shift can bias inferred cosmological pa- 
rameters if not handled carefully. Tests for bias in cosmol- 
ogy constructed by segregating SNe la by host-type were dis- 
cussed in Perlmutter et al. (1997) and carried out with a small 
sample of SNe la (17) in Perlmutter et al. (1999), and a larger 
sample (39) in Sullivan et al. (2003) without any differences in 



fitted cosmological parameters apparent within the statistical 
uncertainties. However, with a much larger sample from the 
SNLS (Sullivan et al. 2010), it becomes possible to see sub- 
tle differences in (3 dependent on the host. SNe la hosted by 
lower specific star-formation rate galaxies obey a shallower 
relation than SNe la hosted by higher specific star-formation 
rate galaxies (/3 ~ 2.8 compared to j3 ~ 3.5; the exact value 
depends on the choice of specific star-formation rate delin- 
eating the low and high star-forming samples). Similarly, us- 
ing SNe drawn from the Sloan Digital Sky Survey (SDSS), 
Lampeitl et al. (2010) find (3 ~ 2.5 for passive hosts and 
f3 ~ 3 for star-forming hosts. Since (3 reflects the combined 
effect of host galaxy dust extinction and an intrinsic SN color- 
luminosity relation, it is not surprising that j3 should be closer 
to the Milky Way dust total-to-selective extinction coefficient 
of Rb = 4.1 in dusty star-forming galaxies than in passive 
galaxies. This result also suggests that as the fraction of low 
star-forming early-type galaxies decreases with redshift and 
the fraction of high star-forming galaxies increases, the aver- 
age value of the relation (f3 (z)) will increase and introduce 
a bias on measurements of w if the subpopulations are not 
treated separately. 

Perhaps a more important consideration, however, is the re- 
cently reported evidence for a correlation between host galaxy 
mass and Hubble residual after applying corrections for light 
curve shape and color. By fitting SEDs derived from the PE- 
GASE2 (Fioc & Rocca-Volmerange 1997, 1999) stellar pop- 
ulation synthesis models to SDSS photometry of host galax- 
ies, Kelly et al. (2010) measured the host masses of the low- 
z SNe in the Hicken et al. (2009) analysis and compared 
these to their Hubble residuals obtained with a variety of light 
curve fitters. They found that SNe la in more massive galax- 
ies are brighter (after stretch and color corrections) by about 
~ 0.1 mag. Sullivan et al. (2010) found a similar relation 
in SNLS SNe and host galaxy photometry using PEGASE2 
SEDs and the SiFTO light curve fitter (Conley et al. 2008). 
Finally, Lampeitl et al. (2010) found a similar relationship in 
the SDSS-II SN survey by segregating SNe by host type, with 
passive galaxies hosting brighter post-correction SNe. How- 
ever, because the passive galaxies in this dataset are on aver- 
age more massive than the star-forming galaxies, this relation 
can also be framed in terms of host galaxy mass. These results 
are particularly important for targeted surveys, such as those 
that target massive galaxies at low redshift or the HST Clus- 
ter SN Survey which targets clusters of massive galaxies at 
high redshift. The demographics of hosts from these surveys 
are substantially different than for untargeted surveys. The 
SNe they discover should be on average brighter (after color 
and light curve shape corrections) than those in other surveys, 
a property which needs to be corrected for in cosmological 
analyses. 

2.4. The HST Cluster SN Survey 

By targeting massive galaxy clusters, the HST Clus- 
ter SN Survey was designed to efficiently discover well- 
characterized SNe la at high redshift. Clusters are rich in 
elliptical galaxies, which constitute a linear red sequence in 
a color-magnitude diagram. The observed evolution of the 
red-sequence slope indicates that ellipticals in the cores of 
clusters have passively evolved since forming at high-redshift 
(z > 2) and that the redder colors of more massive galaxies is 
due to increased metallicity (Kodama 1997; Kodama & Ari- 
moto 1997; Kauffmann & Chariot 1998; Gladders et al. 1998). 
With filters straddling the 4000A break, the red sequence can 
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be readily isolated from foreground and background galaxies, 
a technique also used to find clusters (e.g. Gladders & Yee 
2000). 

Clusters also provide a convenient mechanism for probing 
the dust contents of cluster early-type galaxies. Just as the 
Mg 2 absorption feature traces age and metallicity at low red- 
shift, the magnitude axis of the red sequence traces metallicity 
at high redshift. A dispersion in the residuals from the color- 
magnitude relation of roughly 3% is found in clusters ranging 
from Coma to high redshift (z ~ 1.5) (Bower et al. 1992; Ellis 
et al. 1997; Stanford et al. 1998; Blakeslee et al. 2003a; Mei 
et al. 2006b,a; Lidman et al. 2008; Mei et al. 2009). Some of 
this dispersion can be attributed to differences in galaxy ages, 
and what remains sets a limit on dust. Finally, we also note 
that the uniform old stellar populations of elliptical galaxies 
imply simpler SEDs and hence better mass estimates with 
which to address the trends of §2.3. 

The completed HST Cluster SN Survey produced 19 z > 
0.9 SNe, eight of which were found in the clusters themselves. 
Deep two-color images of the cluster galaxies (including SN 
hosts) were also produced by stacking together SN search and 
follow-up epochs. Similar data exist from the SN surveys of 
the GOODS fields, which have discovered 26 z > 0.9 SNe la, 
all in the field. We analyze the host galaxies of these two sets 
of SNe in this paper. 

3. DATA 

To identify a set of SNe la hosted by a uniform stellar 
population and minimally affected by dust, we compare the 
photometry of SN host galaxies to the photometry of red- 
sequence galaxies and look for star-formation indicators in 
SN host spectroscopy. Identifying the red sequence requires 
careful measurement of galaxy magnitudes and colors, which 
we describe in this section. We also measure quantitative mor- 
phology parameters which we use to enhance the contrast of 
the red sequence in color-magnitude diagrams by exploiting 
the early-type galaxy population dominance along the red se- 
quence. 

3.1. Image reduction and photometry 

Twenty-five massive high-redshift (0.9 < z < 1.46) galaxy 
clusters selected from X-ray, optical, and IR surveys were 
chosen for the HST Cluster SN Survey (Dawson et al. 2009). 
Clusters were each visited by HST four to nine times be- 
tween July 2005 and December 2006. Each visit typically 
consisted of four ~ 500 second exposures in the F850LP 
filter (hereafter zg, 5() ) of the Advanced Camera for Surveys 
(ACS) wide field camera (WFC) and one ~ 500 second ex- 
posure in the F775W filter (hereafter i 775 ) of the ACS WFC. 
The 1775 filter roughly matches rest-frame J7-band for clus- 
ters with 0.9 < z < 1.25, with the best match occuring at 
z = 1.1. The zgso filter roughly matches the rest-frame B- 
band in this redshift range with its closest match occuring at 
z = 1.05. For more distant clusters with 1.25 < z < 1.46, 
the ^850 filter more closely matches rest-frame U -band, with 
the best overlap at z = 1.45. In this paper we look prin- 
cipally at the deep coadditions of exposures from all obser- 
vation epochs. Due to gyroscopic constraints, HST visits to 
individual clusters necessarily occurred at different position 
angles, resulting in coadditions in which pixels near the edges 
have smaller effective exposure time than pixels near the cen- 
ter (Figure 1). Each exposure specifically targeted the cluster 
core, so total integration time is nearly constant in the central 
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FIG. 1.— Final 2 8 50 mosaic for cluster ISCS J1432.4+3332. The differ- 
ences in background noise as a function of position are apparent towards the 
edges of the mosaic. 

region of each coaddition. Four clusters, RDCS J09 10+54 
(Mei et al. 2006a), RDCS J0848+44 (Postman et al. 2005), 
RDCS J1252-29(Blakeslee et al. 2003a) and XMMU 2235.3- 
2557 (Jee et al. 2009), had been previously targeted by ACS 
in i 775 and z S50 (PID9290 and PID9919), and we have in- 
cluded these additional exposures in our coadded images. 
Cluster CL 1604+4304 (Postman et al. 2005) had also been 
observe with ACS, but not in 1775 and zsso- The individually 
sky subtracted exposures were stacked using MultiDrizzle 
(Fruchter & Hook 2002; Koekemoer et al. 2002) with a square 
kernel, pixfrac = 0.8 and the native output pixel scale of 
0.05". 

We similarly processed images obtained by HST ACS as 
part of the HST Great Observatories Origins Deep Survey 
(GOODS) Treasury program (Giavalisco et al. 2004). The 
HST GOODS program targeted two high galactic latitude 
fields to obtain deep multiband images for studies of galaxy 
evolution and included a "piggy-back" SN program to fol- 
low suspected z > 1 SNe with ACS and the Near Infrared 
Camera and Multi -object Spectrograph (NICMOS). The orig- 
inal search consisted of five epochs over 15 ACS pointings for 
each of the two GOODS fields. A subsequent extension con- 
tributed an additional 14 epochs. The survey used nearly iden- 
tical exposures in i 775 and z 850 for SN discovery and follow- 
up as our program. The initial GOODS SN survey yielded 
11 z > 0.9 SNe la and its extension yielded an additional 13 
z > 0.9 SNe la (Riess et al. 2004, 2007). We also add the 
host of one z > 0.9 SN la discovered in a previous search of 
these fields - SN2002dd Blakeslee et al. (2003b). One other 
z > 0.9 SN has been discovered in these fields - SN1997ff, 
but its most likely redshift, z = 1.76 is tentative and far be- 
yond that of the clusters studied here. Lacking a comparison 
cluster red-sequence we do not analyze this SN host. 

Galaxies that hosted SNe required special attention to pre- 
vent the SN light from biasing photometric and morpholog- 
ical statistics of the galaxies themselves. For each SN we 
created postage stamp images of each epoch in both bands. 
From these images, we subtracted a PSF model (described be- 
low) scaled to the flux of the light curve fit (details in Suzuki 
et al. (2011)) of the SN in each epoch. The postage stamp 
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images are then stacked and processed in the same manner as 
the full mosaics. When the final catalogs are created, the val- 
ues derived from the SN-subtracted postage stamp images of 
SN hosts are substituted for the values derived from the full 
mosaics. To test the accuracy of the PSF subtractions, we also 
generated postage-stamp images of SN hosts by stacking only 
the subset of epochs in which the SN was far from its peak. 
Except for a few cases in which only one or two such epochs 
were available, we obtained consistent photometric and mor- 
phological results to those obtained from the PSF-subtracted 
images. 

3.1.1. Object detection 

To form initial object catalogs we used S EXTRACTOR ver- 
sion 2.8.6 (Bertin & Arnouts 1996) in dual image mode, re- 
lying on the 0§5o images for object detection. Determining 
appropriate parameters for the extraction proved challenging 
due to galaxy crowding near cluster cores. One needs to si- 
multaneously be able to deblend neighboring galaxies while 
avoiding the dissection of single complexly structured galax- 
ies into multiple catalog entries. We found that the two-pass 
Cold/Hot method (Rix et al. 2004) produced the best catalogs, 
particularly near cluster cores. This method works by ag- 
gressively deblending relatively bright objects (such as cluster 
core galaxies) in an initial pass of S EXTRACTOR, and subse- 
quently filling in the object catalog with fainter objects in a 
second pass with the deblending parameters set less aggres- 
sively. The specific parameters for the two steps were opti- 
mized by trial and error judged by the successful identification 
and segmentation of galaxies near cluster cores. 

3.1.2. PSF construction 

In order to accurately measure the colors of galaxies we 
require PSF estimates for both the i 77 $ and 2350 bands. To 
identify stars in each field with which to construct a PSF we 
used the SExtractor FLUXJRADIUS statistic. Objects 
with FLUX_RADIUS in the range [1.4, 1.72] pixels are likely 
to be stars. For fields with many stars, we selected all stars 
within 1.3' of the cluster center and zs^q MAG^MJTO be- 
tween 18 and 24. For more sparsely populated fields, we se- 
lected the 10 stars closest to the cluster center falling in the 
same magnitude range. For each band of each field, an initial 
PSF was constructed by subtracting a local background (de- 
scribed below) from each star, oversampling each pixel 9x9 
times with the IDL procedure CONGRID, renormalizing each 
oversampled image by the flux near the stellar core, aligning 
the oversampled images to the nearest subpixel and taking the 
subpixel by subpixel median image. This PSF model was then 
refined by iteratively fitting it to the original stars and adding 
the oversampled median residual back into the model. The fi- 
nal PSF was trimmed to 31 x 31 pixels, which contains more 
than 92% of the encircled energy in both 4775 and 2350 bands 
(Sirianni et al. 2005). 

For more shallowly exposed fields, the wings of the PSF 
derived this way are relatively noisy and tend to be biased high 
when compared with PSFs from more deeply exposed fields. 
To reduce this bias, we match the radial profile of each PSF 
to the radial profile of the PSF derived from a reference field, 
that of cluster RDCS J 1252.9-2927, which was both exposed 
more deeply than other fields and contained many stars. To 
measure the radial profiles of each PSF, we fit a ID B-spline 
to the subpixel values as a function of radius. We then scale 
the PSF subpixel values by the ratio of the target field B-spline 



and the reference field B-spline. This procedure produces PSF 
grids that are consistent with the encircled energy functions of 
Sirianni et al. (2005) to ~ 1%. They also individually retain 
the correct azimuthal structure for coadditions of exposures 
which occured at many different position angles. Finally, we 
note that while a more precise PSF model would vary across 
the field of each exposure, ours is sufficient for correcting the 
effects of differential blurring between the 1775 and 2350 filters 
when computing a galaxy color. 

3.1.3. Galaxy magnitudes and colors 

Nonstellar objects with 19 < z g50 < 26 (measured with 
MAG_AUTO) and -1 < i 775 - z 850 < 2 (measured with 
MAG^VPER with a 10 pixel radius), were then selected as 
galaxies of interest for further processing. For each of these 
galaxies, postage stamp images were cut out of the main ZS50 
mosaic for processing with GALFIT version 3.0 (Peng et al. 
2010). We used GALFIT to obtain two key statistics: the 
Z850 magnitude, and the half-light radius R e . We also used 
GALFIT to produce an interloper-subtracted image used later 
when measuring quantitative morphology parameters. 

Initial parameters for GALFIT were chosen using a variant 
of the GALAPAGOS algorithm (Haussler et al. 2007). Galax- 
ies were modeled as Sersic (1968) profiles with their Sersic 
indices constrained between one and four. Interloper galaxies 
near each target galaxy were either masked out or simultane- 
ously fit depending on the degree of overlap. For each galaxy, 
an adaptively sized elliptical annulus (typically with 6" semi- 
major axis and ~ 3" width) with all galaxies and stars aggres- 
sively masked out was used to find the local sky background 
level, which was held fixed during the fit. We provided GAL- 
FIT with our estimated 2350 PSF to internally convolve with 
its model before fitting to the actual image. The Sersic profile 
magnitudes (and magnitude uncertainties) from the GALFIT 
fits, corrected for Milky Way extinction with the dust maps of 
(Schlegel et al. 1998), are the z S50 magnitudes and uncertain- 
ties reported throughout the rest of this paper. 

Measurements of galaxy color are complicated by the fact 
that the z S50 PSF is ~ 10% broader than the i 775 PSF. To 
account for this difference we implemented a PSF matching 
scheme where, when measuring color, the 2:350 image is con- 
volved with the i 775 PSF and vice-versa. To fairly treat any 
color gradients in galaxies of different angular sizes, we mea- 
sure color within a circular aperture with radius equal to the 
apparent (as opposed to intrinsic) half-light radius R e of the 
cross-convolved 2350 image. Measuring this radius requires 
us to run GALFIT a second time on each galaxy, this time us- 
ing the cross-convolved 2350 image and not providing a PSF 
for GALFIT to internally convolve with its model. A mini- 
mum color aperture radius of 3 pixels, corresponding to 1.2 
kpc at 2 ~ 1, was enforced. 

We modeled the errors for each band's flux as the sum in 
quadrature of the Poisson error from the object flux and the 
error in the background (e.g. sky, CCD readout noise). Be- 
cause some of the clusters are not evenly exposed over the 
analysis region, we estimated the background contribution to 
the error for each galaxy using <7bkg, the local pixel-by-pixel 
standard deviation within the same elliptical annulus previ- 
ously defined in order to measure the local background level. 
If the pixels in the coadded images were uncorrected, then the 
aperture uncertainty er apor turc would be related to the pixel- 
by-pixel uncertainty Obkg by a factor of the square-root of the 
number of pixels in the aperture. However, since both Mul- 
TlDRlZZLE and cross-convolution introduce correlations be- 
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FIG. 2. — Left: Empirically measured photometric error contribution from 
background fluctuations in the 1775 filter, including sky level, readout noise, 
and MULTlDRIZZLE artifacts. The x-axis is the average of o^kg for galaxies 
in the central region of each GOODS tile or cluster field, where a"bkg> as 
described in the text, is the standard deviation of background pixels (in the 
pre - cross-convolved images) surrounding an individual galaxy. The y-axis 
is the standard deviation of circular sky apertures (in the cross-convolved 
images) of fixed radius in the central region of each GOODS tile or cluster 
field. The different curves from bottom to top are computed for aperture radii 
of 4.7, 8.4, and 15.0 pixels; solid red for cluster fields and dashed blue for 
GOODS tiles. The empirical relation is slightly different for the GOODS tiles 
compared to the cluster fields, presumably due to a different mosaic pattern 
and subsequent MULTlDRIZZLE artifacts. Right: The slope of the relation 
in the left-hand panel scaled by 1/ ^JN v i x as a function of aperture radius. 
The solid red curve indicates cluster fields and dashed blue curve indicates 
GOODS tiles. For uncorrected pixels, this relation should be horizontal. The 
deviation from horizontal confirms that the pixels are correlated, justifying 
our empirical calibration. 

tween nearby pixels, we need to empirically calibrate the re- 
lation between cr ap erture and Cbkg- We measured <7bkg (in the 
unconvolved coadded images) and the local background sub- 
tracted flux (in the cross-convolved coadded images) of 1000 
sky apertures (selected to avoid galaxies and stars) in each 
cluster field and GOODS tile. For each aperture center we 
measured the flux at 15 different radii logarithmically spaced 
from 1 pixel to 15 pixels. For a given radius, the average value 
of (Tbkg for each field is tightly correlated with the standard 
deviation of the 1000 aperture fluxes for that field (Figures 2, 
3). The correlation is slightly different between the cluster 
fields and the GOODS tiles, presumably due to different aver- 
age exposure times, different mosaic patterns and subsequent 
MULTlDRIZZLE artifacts. The right-hand panels of Figures 
2 and 3 show the deviations from the scalings expected for 
uncorrelated noise, which would show up as horizontal lines. 

To apply these data to galaxies we interpolate the slope 
of the (Taperture - "bkg relation from the data plotted in the 
right-hand panels of Figures 2 and 3 at the measured galaxy 
half-light radius R e . In cases where R e is larger than 15 pix- 
els, we extrapolate using a linear fit to the data with radius 
between 8 and 15 pixels. This leads to an over-estimate of 
the background uncertainties for very bright galaxies. Fortu- 
nately, less than 6% of the red-sequence galaxies in our sam- 
ple have R e > 15 pixels, and only 1% have R e > 25 pixels. 
In this way the background error is dependent on each individ- 
ual object's exposure depth (quantified by (Tbkg), even when 
the exposure depth varies across the image. Errors in flux are 
then converted to errors in magnitude, and the cataloged color 
error is the sum in quadrature of the 1775 magnitude error and 
the Z850 magnitude error. Exposure depth varies from field to 
field, but typical uncertainties in 1775 — zgso (f° r red galaxies) 
are 0.01 mag at zgso = 21 mag and 0.05 mag at z 850 = 24 
mag. 

3.2. Quantitative morphology 



FIG. 3. — Same as Figure 2 but for the zgso filter. 

We rely on quantitative morphology measurements to iden- 
tify likely red-sequence members. The two parameters we use 
are the Gini coefficient (Abraham et al. 2003) and asymmetry 
(Abraham et al. 1996). The Gini coefficient is a measure of 
the inequality in the distribution of pixel fluxes belonging to 
a galaxy, independent of their positions. The minimum Gini 
coefficient of zero indicates perfect equality; i.e. the flux in 
each pixel is identical. The maximum Gini coefficient of one 
indicates maximum inequality, in which the total flux of all 
pixels within the chosen aperture is actually contained in just 
one pixel (and all other pixels contain zero flux). The con- 
centrated cores of elliptical galaxies typically generate higher 
Gini coefficients for these galaxies than in later-type galaxies. 

The asymmetry measures how poorly a galaxy image 
matches itself when rotated 180 degrees about its center. It 
is defined as half the sum of the absolute values of the 180 
degree subtraction residual pixels divided by the sum of the 
unsubtracted pixels in a given aperture. Elliptical galaxies 
typically have smaller asymmetry than later-type galaxies. 

Just as the crowding in cluster cores is an issue when form- 
ing S Extractor catalogs, it is also an issue when selecting 
apertures for morphology measurements. To overcome this 
difficulty, we construct interloper-subtracted images by using 
GALFIT to subtract models of nearby interfering galaxies for 
each target galaxy. For each interloper-subtracted image, a 
preliminary aperture is created by collecting all pixels exceed- 
ing 1.5fibkg contiguous to the pixel at the target galaxy's cen- 
ter. Using this aperture, we then compute the galaxy's quasi- 
Petrosian flux (Abraham et al. 2007). If the calculation con- 
verges, then a new isophotal aperture is created by selecting 
pixels from the preliminary aperture which exceed the quasi- 
Petrosian flux. The quasi-Petrosian flux fails to converge for 
a small number of galaxies which are then excluded from the 
subsequent analysis. These galaxies are usually quite faint 
and would have had large uncertainties had the fits converged. 

The Gini coefficient is measured within the isophotal aper- 
ture. Since the Gini coefficient calculation does not depend 
on the relative geometry of the aperture pixels but only on the 
pixel values, the error can be estimated from bootstrap resam- 
pling (Abraham et al. 2003) (i.e. by repeatedly resampling 
from the original pixel values, with replacement, and recom- 
puting the Gini coefficient on these resamplings. We recom- 
pute the Gini coefficient this way 1000 times to determine 
the Gini coefficient probability distribution function for each 
galaxy and record the standard deviation of this distribution 
as the Gini coefficient error). 

We measure the asymmetry within a symmetrized aper- 
ture consisting of the intersection of the isophotal aperture 
with its 180 degree rotation. The center of the rotation is fit 
iteratively to minimize the asymmetry (so the symmetrized 
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aperture changes with each iteration, depending on the cur- 
rent candidate center). Since the random fluctuations present 
in the background contribute some asymmetry (positive by 
definition), a correction must be applied. To estimate this 
correction, we generate 1000 Gaussian background images 
with standard deviation Obkg an d measure their asymmetry 
(but without normalization) in the same symmetrized aper- 
ture as the galaxy. We can then subtract this background con- 
tribution to the galaxy's asymmetry. The standard deviation 
of these background measurements is our estimated error for 
each galaxy's asymmetry (note that galaxy photon noise is 
much smaller than the background noise). 

3.3. Spectroscopy 

As the HST Cluster SN Survey produced SN candidates, 
they were spectroscopically targeted using prescheduled ob- 
serving time on DEIMOS on Keck II (Faber et al. 2003) and 
FOCAS on Subaru (Kashikawa et al. 2002) and with ToO 
requests on FORS1 and FORS2 on Kueyen and Antu at the 
VLT (Appenzeller et al. 1998). The FORS1, FORS2 and 
DEIMOS observations are described in Dawson et al. (2009); 
the FOCAS observations are described in Morokuma et al. 
(2010). By observing SN candidates through slitmasks we 
were able to simultaneously target likely cluster members and 
form spectroscopic catalogs of many galaxies in the cluster 
fields. Galaxy redshifts were found through cross-correlation 
with template eigenspectra derived from SDSS spectra (Agol 
etal. 2011). 

The flux and equivalent width of the [O II] 3727A emission 
line doublet was measured by fitting the simple stellar popula- 
tion (SSP) templates from Bruzual & Chariot (2003, herefac- 
ter BC03) assuming a Chabrier (2003) initial mass function to 
data on both sides of the feature to define the continuum. The 
flux and equivalent width were computed by integrating from 
500 km/s blueward of the line centered at 3726. 032A to 500 
km/s redward of the line centered at 3728. 815A. To mitigate 
slit losses in the spectroscopy, we normalized the spectra to 
the observed 1775 photometric magnitude. 

Additional literature redshifts and [O II] equivalent width 
measurements were used as available (Andreon et al. 2008; 
Bremer et al. 2006; Brodwin et al. 2006; Demarco et al. 2007; 
Eisenhardt et al. 2008; Hilton et al. 2007, 2009; Postman et al. 
1998; Rosati et al. 1999; Stanford et al. 2002, 2005). Spectra 
of z > 0.9 SN hosts from the HST Cluster SN Survey are 
presented in Figure 4. A few of these spectra merit individual 
discussion. The redshift of SN SCP06T1 is determined from 
a single well-detected emission line, which we assume to be 
[O II]. Likewise, the redshift of SN SCP06X26 is determined 
from a single [O II] emission line, however the detection of 
this line is much more tentative than for SN SCP06T1 and 
is only barely visible in the 2D spectrum. Finally, we have 
not identified any features in the spectrum of SN SCP06E12 
that would allow us to determine a redshift. The color of 
this galaxy is consistent with the color of the red sequence 
of the targeted cluster in the same field of view, so it may be 
a cluster member. However, we note that the colors of the 
four early-type field SN hosts from the HST Cluster SN Sur- 
vey identified later are also roughly consistent with the colors 
of the red sequences of the targeted clusters in their fields of 
view. Though there is considerable uncertainty in the redshift 
of SN SCP06E12 and its host, we follow through with our 
analysis of this galaxy assuming it is a cluster member. We 
note that it is not used for cosmological analyses or analyses 
of SN correlations with their hosts. All of the remaining SN 



hosts have secure redshifts. 

4. ANALYSIS 

In this section we describe our determination of the red 
sequence in clusters and in the field, which we later use to 
characterize SN host galaxies. To consistently identify cluster 
red-sequence galaxies, especially in poorer clusters where the 
contrast between bluer foreground or late-type galaxies and 
the cluster red sequence is less distinct, we generate a com- 
posite red sequence using spectroscopically confirmed mem- 
bers of multiple clusters with colors loosely consistent with 
the red sequence. When there is no cluster red sequence with 
which to compare SN host galaxies, as is the case for the 
field hosted SNe in the HST Cluster SN Survey and for all the 
GOODS fields SNe, this composite red sequence will serve as 
the reference red sequence. We use an evolving SED model to 
^-correct and evolution correct red-sequence members from 
their original redshifts to a target redshift. In this way we pre- 
dict the location of the red sequence across the redshift range 
of the clusters in the survey. 

4.1. Quantitative morphology calibration 

Cluster red sequences are dominated by early-type galax- 
ies. To enhance the contrast of the red sequence against 
the background color magnitude diagram, we use the 
Gini coefficient and asymmetry from §3.2 to select these 
galaxies. We calibrate our selection against visually de- 
rived morphology catalogs (by co-author MP) for three 
clusters: ISCS J1432.4+3332, RDCS 1252.9-2927, and 
ISCS J1438.1+3414 at redshifts 1.10, 1.24, and 1.41 respec- 
tively. Figure 5 shows the locus of red (1775 — 2:350 > 0.75) 
visually classified elliptical galaxies compared to other visu- 
ally classified morphological types in the Gini coefficient - 
asymmetry plane. The color cut eliminates many of the bluer 
galaxies (which we are not interested in here) but is not so 
red (for these clusters) that red-sequence members are cut. To 
select visually classified ellipticals we can evidently choose 
galaxies with Gini coefficient greater than 0.45 and asymme- 
try less than 0.08, cuts which we refer to as our narrower 
morphological cuts. With these cuts, ~ 90% of all red vi- 
sually classified elliptical galaxies are selected. SO galaxies 
make up about 40% of the objects selected. The contamina- 
tion from late-type galaxies is < 5%. Later, when we fit the 
CMRs of clusters we choose to use somewhat more lenient 
cuts of Gini > 0.4 and asymmetry < 0.1, which we refer to 
as our broader morphological cuts. This increases the num- 
ber of red-sequence early-type galaxies included by ~ 50% 
but only increases the late-type contamination to ~ 10%. The 
broader cuts further increase the contrast of the red sequence 
but, critically, we find they only marginally increase the scat- 
ter of the red sequence (from 0.045 ± 0.003 to 0.048 ± 0.003 
when clusters are stacked; §6). 

As an independent test of the reliability of our cuts, we also 
investigate the strength of [O II] line emission in the Gini co- 
efficient - asymmetry plane. Figure 6 shows the locus of red 
galaxies with weak or no [O II] emission (equivalent width 
greater than — 5A). Our broader cuts select ~ 90% of [O II] 
quiet galaxies, though the contamination of [O II] emitting 
galaxies is higher (~ 35%) than the contamination of late- 
type galaxies in the visual morphology case. This is partly a 
selection effect: redshifts are easier to obtain from emission 
lines than from absorption lines so [O II] emitting galaxies 
are more likely to be spectroscopically confirmed as cluster 
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FIG. 4. — SCP supernovae host spectroscopy. For each host, data are in black, and the best fitting linear combination of eigenspectra at the fitted redshift is in 
green. The blue arrow indicates the wavelength of the [O II] 3727 emission line doublet. Note that the green line, used here for determining the redshift, is not 
the same as the BC03 template (which does not contain [O II] emission) used in fitting the [O II] equivalent width. 
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FIG. 5. — Comparison of visually derived morphologies to quantitative 
morphology measurements for three clusters. Blue boxes represent visually 
classified late-type galaxies with T-type from to 6. Green triangles repre- 
sent SO galaxies with visual T-type from —2 to —1. Red circles represent 
elliptical galaxies with visual T-type from —6 to —3. Colored filled symbols 
represent galaxies with 2775 — 2850 > 0.75, consistent with the red sequence, 
light grey open symbols represent galaxies bluer than the red sequence. The 
two boxes in the lower-right show the position of our broader and narrower 
morphology cuts. The horizontal and vertical lines extending from plotting 
symbols indicate lcr uncertainties. 

members and hence appear in the plot. The observed fraction 
of [O II] emitting galaxies is compatible with results reported 
in Postman et al. (2001), in which 45% of galaxies within 
the central 1.0/igg 1 Mpc were observed to have [O II] equiv- 
alent width less than — 15A. Also, Yan et al. (2006) have an- 
alyzed ~ 55000 SDSS low-redshift spectra and determined 
that while nearly 38% of red galaxies show [O II] emission, 
only 9% of these show the emission line ratios characteris- 
tic of star formation. The remaining [O II] -emitting galaxies 
show line ratios indicative of LINERS or AGN. Lemaux et al. 
(2010) have recently extended this analysis to high-redshift 
by obtaining near-infrared spectroscopy of z > 0.8 cluster 
galaxies. Of the five [O ll]-emitting red-sequence early-type 
galaxies in their sample, four were found to be consistent with 
AGN/LINERS, with some evidence that the fifth might con- 
tain AGN/LINER activity as well. We are unable to directly 
check for AGN/LINER-like line ratios in our spectra because 
at high redshift the necessary comparison lines are shifted out 
of the wavelength range of our spectroscopic coverage. How- 
ever, we suspect that most of the [O II] -emitting morpholog- 
ically early-type red galaxies in our sample are not actively 
forming stars. 

4.2. The composite red sequence 

The identification of red-sequence galaxies in clusters with 
incomplete spectroscopy can be challenging, especially for 
poorer clusters where it is difficult to distinguish red-sequence 
members from interloping galaxies in the foreground or back- 
ground of the cluster. We therefore use the spectroscopically 
confirmed red-sequence members from all of our clusters to 
construct a composite red sequence which serves as both a 
guide for fitting the CMR of individual cluster red sequences 
and also as a reference red sequence for field elliptical SN 



FIG. 6. — Spectroscopically confirmed galaxies' asymmetry and Gini co- 
efficient distributions. Red circles indicate galaxies with an [O II] equiva- 
lent width greater than — 5A (small or no emission). Blue squares indicate 
galaxies with strong [O II] emission with an equivalent width less than — 5A. 
Colored filled symbols represent galaxies with 4775 — 2850 > 0.75, consis- 
tent with the red sequence, light grey open symbols represent galaxies bluer 
than the red sequence. The two boxes in the lower-right show the position of 
our broader and narrower morphology cuts. The horizontal and vertical lines 
extending from plotting symbols indicate lcr uncertainties. 

hosts. We start by assuming that our sample of spectroscop- 
ically confirmed cluster members with quantitative morphol- 
ogy measurements consistent with elliptical galaxies (passing 
our broader morphology cuts) and colors loosely consistent 
with an old stellar population are each members of their own 
cluster's red sequence. By using a suitable ^-correction and 
evolution correction, we can then project each galaxy's mag- 
nitude and color from their measured values at the galaxy's 
original redshift to a target redshift. To accomplish this we 
assume that the each sample galaxy's SED takes the form: 

A/ gal BC03 A (T(z),Z,z) 
M Q {D L {z)f 

where Dl(z) is the luminosity distance to redshift z, 
BC03a (T (z) , Z, z) indicates the luminosity density of the 
BC03 SSP template SED of age T (z) and metallicity Z with 
wavelength redshifted by (1 + z), and M ga i is the initial stel- 
lar mass of the BC03 template. We only have two band pho- 
tometry for galaxies, so we can only constrain two parame- 
ters affecting the SED. One parameter needs to be the size 
of the galaxy (M ga i), leaving the choice of either constrain- 
ing the galaxy age or metallicity as the second parameter. 
Since these two parameters are largely degenerate, we sim- 
ply fix the age T(z) such that galaxies all form at the same 
redshift Zf OTm and fit for the metallicity Z. In a few cases, 
the range of colors produced by Equation 1 for the available 
BC03 metallicities was not red enough to cover the observed 
galaxy color. This could happen either because of a chance 
positive fluctuation of the galaxy color or because the partic- 
ular galaxy happens to have formed earlier than our chosen 
Zfoim- Since we are not interested in the metallicity values 
themselves, but only in using them to project colors and mag- 
nitudes from one redshift to another, we allow this parameter 
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FIG. 7. — The composite red sequence at redshifts 1.0 (blue diamonds), 
1.2 (green circles), and 1.4 (red squares), with linear fits. The size of the 
plotting symbol is proportional to the weight of the galaxy in the linear fit. 
The dashed lines indicate the measured residual scatter about the fits. The 
black line shows the predicted evolution from z = 0.9 to z = 1.5 of an L* 
cluster member. The redward then blueward evolution of the red sequence 
with redshift is the result of the rest-frame 4000A Balmer break crossing the 
gap between the 4775 and z$rM filters. 

to be extrapolated outside of the BC03 range. In practice, 
we extrapolated metallicities for 27 out of 198 galaxies to 
a maximum of Z = 0.07 (the range of BC03 metallicities is 
0.0001 to 0.05) when choosing Zf orm = 3.0. To then project 
each galaxy to a target redshift, we evaluate the best-fit BC03 
template at age T (z) given the target redshift z, which au- 
tomatically captures luminosity evolution. We calculate 1775 
and zg5o from this template. This procedure is imperfect, and 
in particular can introduce errors when ^-correcting and evo- 
lution correcting over large differences in redshift. To miti- 
gate this effect we de-weight galaxies with redshifts far from 
the target redshift. Specifically, we assign a Gaussian weight 
w oc exp(— (z — Ztarg) 2 /(2 • 0.15 2 ) to each galaxy. The 
characteristic size of the redshift window in this weighting, 
0.15, is arbitrary but we have confirmed that our results do 
not significantly change with window sizes ranging from 0. 10 
to 0.25. The composite red sequences at redshifts 1.0, 1.2, 
and 1.4 are shown in Figure 7 for Zf OTm = 3.0. 

4.3. Individual cluster red sequences 

To determine the red sequences of individual clusters we as- 
sume a model in which galaxies are normally distributed with 
intrinsic scatter <7i nt in color about the mean color-magnitude 
relation (characterized by slope m and intercept at z 850 = 22, 
622) and have their magnitudes distributed like a Schechter 
(1976) function (described by parameters </>*, L* and a). Ad- 
ditionally, foreground and background galaxies will contam- 
inate the color-magnitude diagram. We model the density of 
these galaxies in color-magnitude space as a plane character- 
ized by its average value (B) and derivatives with respect to 
magnitude and color (dB/dzs5o and dS>/d(in<--, — ^sao))- To 
break degeneracies between the red-sequence parameters and 
the foreground/background parameters, we simultaneously fit 
each cluster field with the GOODS fields which serve as a 
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FIG. 8. — Comparison of the slopes of the composite red-sequence CMR 
fits to individual cluster CMR fits when allowing the slope to float. The multi- 
color curves are fits to the composite red sequence assuming different values 
of Zf orm in the range [2.5,7.0]. The black points show the locations of maxi- 
mum likelihood from the MCMC fits to 23 individual clusters. The errorbars 
are the smallest confidence intervals containing 68% of the marginal posterior 
probability. 



control. 

We consider only galaxies that pass our broader mor- 
phology cuts and fall within a redshift dependent color- 
magnitude rectangle for this analysis. The allowed color 
range is &22 mp ±0.3 mag where 622™ P i s the intercept of 
the composite red sequence. The allowed magnitude range 
^ [M* odol - 2.2, M* odel + 1.0] where M^ odcl is determined 
by passive evolution (Eqn. 1) of a red-sequence characteris- 
tic magnitude of 2350 = 22.7 at redshift z = 1.24 (Blakeslee 
et al. 2003a). We construct a likelihood function for these 
galaxies and the GOODS fields control galaxies falling in the 
same color-magnitude rectangle as in Andreon et al. (2005) 
and use a Markov chain Monte Carlo algorithm to explore the 
nine-dimensional space (<7j n t, m, 622, 4>* \ L*, a, B, dB/dzsso 
and dB/d(i775 — Zsso))- The techniques presented in Dunk- 
ley et al. (2005) are used to determine the convergence of the 
Markov chains. The marginal distributions of fitting parame- 
ters can then be extracted directly from the Markov chains. 

In Figure 8, we compare the slopes of individual cluster red- 
sequence CMRs to the slopes of the composite red-sequence 
CMRs generated assuming a variety of values of Zf oim be- 
tween 2.5 and 7.0. The effects on the composite CMR slope 
of changing Zf orm are much smaller than the statistical uncer- 
tainties in individual CMR slope fits. In fact, the slopes of 
the individual cluster fits are statistically consistent with the 
narrow range of slopes of the composite red sequence regard- 
less of Zf orm (\ 2 = 19 for 23 clusters assuming Zf orm = 3 
and changes negligibly for other values of Zf onl ). In order to 
consistently treat both field and cluster hosted SNe and also to 
avoid using physically perplexing positive slopes, we choose 
to fix the slope of each cluster to the value of the composite 
red sequence. We choose Zf olm = 3.0 as a typical galaxy for- 
mation epoch and refit each cluster red-sequence CMR. The 
results of these fixed-slope fits are presented in Table 1 . 
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TABLE 1 
Color-Magnitude Relation fits 



ID 



Cluster Name 



Redshift Intercept Slope 0(1775 — zsso) c(£/ — V) z 



al 



A 


XMMXCS J2215.9-1738 


1.457 


0.923 


-0 


.020 





.051; 


B 


XMMU J2205.8-0159 


1.12 


1.028 


-0. 


.023 





.040; 


C 


XMMU J1229.4+0151 


0.974 


0.906 


-0 


.028 





.066; 


D 


RCS J022144-0321.7 


1.017 


0.953 


-0 


.026 





.069; 


E 


WARPS J1415. 1+3612 


1.026 


0.969 


-0 


.025 





.060; 


F 


ISCS J1432.4+3332 


1.110 


1.031 


-0. 


.023 





.046; 


G 


ISCS J1429.3+3437 


1.259 


0.989 


-0 


.019 





.062; 


H 


ISCS J1434.4+3426 


1.241 


0.964 


-0 


.018 





.061 


I 


ISCS J1432.6+3436 


1.347 


0.969 


-0 


.023 





.042; 


J 


ISCS J1434.7+3519 


1.37 












K 


ISCS J1438. 1+3414 


1.414 


0.944 


-0 


.021 





.008; 


L 


ISCS J1433. 8+3325 


1.37 












M 


CL J 1604+4304 


0.897 


0.749 


-0 


.026 





.050; 


N 


RCS J022056-0333.4 


1.026 


0.982 


-0. 


.025 





.043; 


P 


RCS J033750-2844.8 


1.1" 


0.988 


-0 


.023 





.023; 


Q 


RCS J043934-2904.7 


0.955 


0.858 


-0 


.027 





.045; 


R 


XLSS J0223.0-0436 


1.215 


1.013 


-0 


.019 





.015 


S 


RCS J215641-0448.1 


1.07 


1.022 


-0. 


.024 





.005; 


T 


RCS2J151104+0903.3 


0.971 


0.899 


-0 


.028 





.087; 


U 


RCS J234526-3632.6 


1.037 


0.957 


-0 


.025 





.047' 


V 


RCS J23 1953+0038.0 


0.900 


0.766 


-0 


.026 





.051 


w 


RX J0848.9+4452 


1.261 


0.985 


-0. 


.019 





.028; 


X 


RDCS J0910+5422 


1.101 


1.001 


-0 


.023 





.060; 


Y 


RDCS J1252.9-2927 


1.237 


0.965 


-0 


.018 





.029; 


Z 


XMMU J2235.3-2557 


1.390 


0.926 


-0 


.022 





.031; 



+0.034 
-0.020 
+0.186 
-0.002 
+0.003 
-0.014 
+0.017 
-0.017 
+0.015 
-0.013 
+0.016 
-0.016 
+0.019 
-0.034 
+0.037 
-0.013 
+0.014 
-0.025 



+0.011 
-0.008 



+ 0.01!) 
-0.006 
+0.014 
-0.004 
+0.026 
-0.010 
+0.021 
-0.010 



+0.084 
-0.005 



-0.007 
+0.027 
-0.007 
+0.013 
-0.005 
+0.024 
-0.012 
+0.003 
-0.007 
+0.017 
-0.016 



0.101 
0.079 
0.131 
0.137 
0.119 
0.093 
0.125 
0.122 



0.016 



0.086 



0.174 
0.093 



+0 


008 


-0 


.040 


1+0 


372 


'-0 


.003 


+0 


007 


-0 


.028 


•+0 


034 


-0 


.035 


i+0- 


020 


'-0 


.026 


,+0 


032 


'-0 


.031 


+0 


.038 


'-0 


.068 


,+0 


073 


-0 


.026 


,+0 


020 


'-0 


.050 


;+° 


023 


'-0 


.016 


i+° 


.038 


-0 


.012 


;+° 


027 


'-0 


.008 


■+o 


.051 


-0 


.020 


+0 


.041 


-0 


.010 




034 


<t°o 


.020 




.168 


'±°o 


.010 


+0 


080 




.045 


■+o 


030 


'-0 


.014 


+0 


.054 




.015 


■+o 


026 




.000 


l+° 


047 


-0 


.024 


.+0 


007 


'-0 


.014 


+0 


.033 




.032 



23.3 
9.7 
54.8 
32.1 
27.0 
18.8 
23.3 
22.5 
14.8 

30.4 

27.1 
30.4 
15.6 
24.8 
20.6 
7.1 
15.7 
29.1 
30.3 
20.1 
18.4 
34.5 
17.4 



NOTE. — Intercept indicates the CMR color at 2:350 — 22. The slope for each fit is fixed to the slope of the composite 
red sequence at the cluster redshift. (7(2775 — z 85o) (^(U — V) z= q) is the observer-frame (rest-frame) intrinsic scatter of 
morphologically selected red-sequence members. The uncertainties reported indicate the smallest intervals containing 68% of the 
posterior probability for £7(1775 — 2350) and a(U — V) z =q. Af ga l is the effective number of red-sequence galaxies in the fitted 
magnitude range computed by integrating a Schechter function with the best fit values of 0*, L* , and a over the magnitude range 
of the fit. 

a Photometric cluster redshift. 



5. SN HOST CLASSIFICATION 

In this section, we describe how, for each SN host, we use 
the color-magnitude, quantitative morphology, and spectro- 
scopic information to classify the host as either a passively 
evolving early-type galaxy, or a potentially star-forming late- 
type galaxy. We restrict our attention to redshifts greater than 
0.9, which is the lowest cluster redshift in the HST Cluster 
SN Survey. Dawson et al. (2009) lists 17 SNe with z > 0.9. 
From this list we subtract SN SCP06C1, a spectroscopically 
confirmed SN la at z — 0.98 with an uncertain host, but for 
completeness add SN SCP06E12 and SN SCP06X26 for con- 
sideration, which are likely z > 0.9 SNe but lack conclusive 
spectroscopic redshifts. The host of SN SCP06E12 is faint 
but has photometry consistent with the red sequence of the 
cluster in the same field of view, and we assume here that its 
redshift is that of this cluster: 1.026. Spectroscopy of the host 
of SN SCP06X26 shows a possible emission line at 9100A, 
which if [O II] indicates z = 1.44. There are 26 z > 0.9 
SNe la discovered in surveys of the GOODS fields, 19 of 
which are spectroscopically confirmed as Type la (Gilliland 
et al. 1999; Blakeslee et al. 2003b; Riess et al. 2004, 2007). 



We were unable to identify the host galaxy of one of these: 
HST04Gre. This SN is approximately 1.7" from each of two 
potential hosts. Also, the host of SN1997ff has a tentative red- 
shift much greater than that of any of the clusters considered 
here, preventing us from comparing its color and magnitude 
to the red sequence. We do not consider either of these SNe 
in the following analysis. 

The principle criteria that our hosts must pass to be con- 
sidered passively evolving early-types are: (A) they pass our 
broader morphology cuts and (B) they lie within 2er of the red 
sequence where a = \J 'of nt + er; ol , and er co i is the color un- 
certainty. For (B) we use the cluster red sequence for cluster 
member hosts, and the composite red sequence with an as- 
sumed intrinsic scatter of 0.10 mag (greater than the measured 
intrinsic scatter of any of our clusters) for field hosts. Where 
spectroscopy is available, we also take note of possible star 
formation indicators. 

In summary, from the HST Cluster SN Survey we classify 
six out of nine cluster SN hosts and four out of nine field SN 
hosts as early-type. One host galaxy, that of SN SCP06E12, 
we were unable to classify due to its faintness. The remaining 
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HST Cluster SN Survey hosts are classified as late-type. From 
the GOODS SN Survey, we classify seven out of 24 field SN 
hosts as early-type, and the remaining 17 hosts as late-type. 
Classifications are presented in Tables 2 and 3 and details are 
below. 

Early-type Hosts of SCP SNe 

Cluster members 

SNSCP06H5, SNSCP06K18, SN SCP06R12. The hosts of 
these three SNe la are spectroscopically consistent with early- 
type galaxies, have photometry placing them on their clusters' 
red sequences, and pass our narrower morphology cuts. 

SN SCP05D0. The host of this SN has photometry placing 
it on its cluster's red sequence and passes our broader mor- 
phology cuts. The somewhat large asymmetry of this galaxy 
appears to result from a slightly skewed core light distribu- 
tion, though not from any spiral structure or clumps of star 
formation. A small amount of [O II] emission is observed in 
its spectrum which might indicate star formation but is also 
not uncommon in passively evolving galaxies with LINERs 
in their cores. 

SN SCP06K0. The host of this SN is spectroscopically con- 
sistent with an early-type galaxy, has photometry placing it on 
its cluster's red sequence, and passes our broader morphology 
cuts, though it has a somewhat small Gini coefficient. We 
note that the Gini coefficient is biased small for low signal- 
to-noise ratio images (Lotz et al. 2004) (such as the image of 
SN SCP06K0's host) and it is likely that with increased expo- 
sure time the Gini coefficient would increase. 

SN SCP06U4. The host of this SN has photometry placing 
it on its cluster's red sequence and passes our narrower mor- 
phology cuts. Its spectrum exhibits moderate [O II] emission 
but otherwise is consistent with an early-type spectrum. 

Field galaxies 

SN SCP06C0, SN SCP06D6. The hosts of these SNe have 
photometry placing them on the composite red sequence and 
pass our narrower morphology cuts. Small amounts of [O II] 
emission are observed in their spectra. 

SN SCP06G4. The host of this SN has photometry plac- 
ing it on the composite red sequence and passes our narrower 
morphology cuts. Its spectrum shows relatively strong hy- 
drogen Balmer absorption features consistent with an "E+A" 
spectrum (Dressier & Gunn 1983, 1992), which may indi- 
cate it is younger than the other early-type galaxies in our 
sample (though still likely older than the cutoff timescale for 
SNe CC). 

SN SCP06A4. The host of this SN has photometry plac- 
ing it on the composite red sequence and passes our narrower 
morphology cuts. A small amount of [O II] emission may be 
present in its spectrum and hydrogen Balmer absorption lines 
indicate an "E+A" classification. 

Late-type Hosts of SCP SNe 

Cluster members 

SN SCP06B4. The host of this SN fails our broader mor- 
phology cuts, is significantly bluer than its cluster's red se- 
quence, and shows strong [O II] emission in its spectrum. It 
is likely that this galaxy is star-forming. 



SN SCP06F12. Although the host of this SN only barely 
fails our broader Gini coefficient cut, its color is significantly 
bluer than its cluster's red sequence, and it shows strong [O II] 
emission in its spectrum. 

Field galaxies 

SN SCP05P1, SN SCP06T1. The hosts of these SNe fail 
both of our broader morphology cuts, are significantly bluer 
than the composite red sequence, and show strong [O II] emis- 
sion in their spectra. 

SN SCP06N33. Although the host of this SN passes our 
broader Gini coefficient cut, its color is significantly bluer 
than the composite red sequence and it shows strong [O II] 
emission in its spectrum. 

SN SCP06X26. Although we do not detect strong [O II] 
emission from the host of this SN, it fails our broader Gini 
coefficient cut and its color is significantly blue. 

SN SCP06G3. Although the host of this SN passes our nar- 
row Gini coefficient cut, passes (barely) our broader asym- 
metry cut and is only moderately bluer than the composite 
red sequence, its spectrum shows very strong [O II] emission. 
We also note that the axis-ratio of the zg^o GALFIT model 
for this galaxy (b/a — 0.25) is 3.5c smaller than the model 
axis-ratios of the galaxies classified as early-types above; it is 
likely an edge-on disk galaxy. 

Unclassified Hosts of SCP SNe 

SN SCP06E12. We do not have a redshift for the host of 
this SN. Its color is consistent with the red sequence of the 
cluster in the same field of view. Its Gini coefficient is small 
but this galaxy, like SN SCP06K0, has a low signal-to-noise 
ratio image and its Gini coefficient is likely biased low. We 
do not have enough information to classify this galaxy, and do 
not include it in any subsequent analysis in this paper. 

Host Galaxies of GOODS SNe 

From the GOODS SN Survey, we analyze the morphol- 
ogy and photometry (but not spectroscopy) of 24 of the 26 
z > 0.9 SN hosts mentioned in Riess et al. (2007). The hosts 
of seven SNe pass our morphological and photometric cuts: 
SN2003az, SN2003XX, SN2003es, HST04Sas, HST04Tha, 
HST05Lan and SN2002hp. Riess et al. (2007) classifies the 
hosts of all of these but SN2003az as elliptical. The remain- 
ing 17 GOODS SN la hosts analyzed here are classified as 
late-type. 

6. SN TYPING BY HOST GALAXY 

The classification of the supernovae discovered in the HST 
Cluster SN Survey uses several approaches simultaneously 
(see Barbary et al. (2010) for a discussion of the photomet- 
ric and light curve constraints). Here we discuss the typing 
constraints that come from just the host galaxy information. 
These constraints are already quite strong. 

As described in §2.1, classifying a SN as Type la when 
hosted by an early-type galaxy is a robust alternative to spec- 
troscopic typing. However, it is natural to wonder if misclas- 
sification of the host galaxy and subsequently the SN is more 
frequent at high redshift, where morphological classification 
is more difficult and star formation more prevalent. To in- 
vestigate this potential shortcoming, we estimate the a pri- 
ori relative rates of detecting SNe la and SNe CC in the host 
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TABLE 2 
Early-type z > 0.9 SN Hosts 



Name 


^host ' 


'cluster 


[O II] EW (A) 


Confidence Interval 


Gini Coefficient 


Asymmetry 


7."7TK 7(cA 

L i I O ^OOU 


CMR residual 


Notes 


SCP Cluster SN Hosts 






















SN SCP05D0 


1.014 


1.017 


-1.7 


[ -2.6, 


-0.6] 


0.525 ±0.021 


0.084 ± 0.003 


0.939 ±0.021 


—0.016 ± 0. 


072 


s,h,x,c 


SN SCP06H5 


1.231 


1.241 


-0.4 


[ -1.5, 


—0.0] 


0.541 ±0.017 


0.059 ± 0.003 


1.075 ± 0.022 


±0.106 ± 0. 


065 


L,h 


SN SCP06K0 


1.416 


1.414 


-0.4 


[ —2.4, 


-0.0] 


0.419 ±0.018 


0.027 ± 0.008 


0.951 ± 0.066 


±0.020 ± 0. 


066 


L,h,x,c 


SN SCP06K18 


1.412 


1.414 


-0.1 


[ —0.7, 


—0.0] 


0.456 ±0.015 


0.021 ± 0.006 


0.986 ±0.037 


±0.038 ± 0. 


038 




SN SCP06R12 


1.212 


1.215 


-0.2 


[ -1-1, 


—0.0] 


0.455 ± 0.025 


0.054 ± 0.011 


0.956 ±0.050 


—0.019 ± 0. 


052 


L.h 


SN SCP06U4 


1.050 


1.037 


-9.6 


[ —10.9, 


-8.1] 


0.481 ± 0.018 


0.055 ± 0.004 


0.926 ±0.017 


— O.Uoo ± U. 


050 


s,c 


SCP Field SN Hosts 






















SN SCP06A4 


1.193 


1.457 


-1.6 


[ -3.6, 


-0.0] 


0.468 ±0.021 


0.020 ± 0.006 


0.976 ±0.044 


-0.013 ±0. 


109 


L,h,x,c 


SN SCP06C0 


1.092 


0.974 


-3.6 


[ -5.6, 


-1.7] 


0.476 ±0.021 


0.027 ± 0.003 


1.072 ±0.017 


±0.014 ±0. 


101 


L,h,x,c 


SN SCP05D6 


1.314 


1.017 


-3.2 


[ -4.2, 


-2.2] 


0.516 ±0.017 


0.043 ± 0.005 


0.960 ±0.027 


±0.044 ±0. 


103 


L,h,x,c 


SN SCP06G4 


1.350 


1.259 


-0.3 


[-1.5, 


-0.0] 


0.492 ±0.022 


0.000 ± 0.006 


1.004 ± 0.040 


±0.103 ±0. 


108 


s,h,x,c 


GOODS SN Hosts 






















SN2002hp 


1.305 










0.501 ± 0.025 


0.032 ± 0.005 


0.919 ±0.019 


±0.015 ±0. 


102 


h,x,c 


SN2003az 


1.270 










0.435 ± 0.023 


0.054 ± 0.008 


0.927 ± 0.048 


±0.018 ±0. 


111 


s,h,x,c 


SN2003es 


0.954 










0.497 ±0.018 


0.047 ± 0.001 


0.942 ± 0.006 


±0.039 ± 0. 


100 


s,x,c 


SN2003XX 


0.935 










0.504 ±0.028 


0.025 ±0.001 


0.789 ± 0.005 


-0.081 ± 0. 


100 


s,h,x,c 


HST04Sas 


1.390 










0.431 ± 0.023 


0.071 ± 0.007 


0.792 ± 0.032 


-0.081 ± 0. 


105 


s,h,x,c 


HST04Tha 


0.954 










0.432 ± 0.031 


0.058 ±0.001 


0.886 ± 0.008 


±0.027 ±0. 


100 


s,x,c 


HST05Lan 


1.235 










0.492 ± 0.022 


0.078 ± 0.004 


0.997 ± 0.016 


±0.050 ± 0. 


101 


s,h,x,c 



NOTE. — The [O II] EW reported is the median of the posterior probability of the true EW given the observed EW, the spectroscopic uncertainties, and a prior that the EW be positive. 
The reported confidence interval is the smallest interval containing 90% of the posterior probability. The CMR residual is computed using the fixed-slope fits to individual clusters for 

cluster member hosts and the composite red sequence fits for field hosts. The CMR residual uncertainties reported are J cr? nt + cr^ ol where (7i n t is the measured intrinsic scatter of the 
parent cluster for cluster member hosts and a conservative 0.10 mag for field hosts. 

8 SN is classified as Type la from spectrum (Blakeslee et al. 2003a; Riess et al. 2004, 2007; Barbary et al. 2010) 
L SN is classified as Type la from light curve (Riess et al. 2004, 2007; Barbary et al. 2010) 
^ SN passes Union2. 1 cuts and has reliable Hubble residual measurement (Suzuki et al. 20 1 1 ) 
x SN has SALT2 XI uncertainty less than 1.0 (Suzuki et al. 201 1) 
c SN has SALT2 color uncertainty less than 0.1 (Suzuki et al. 201 1) 



galaxies from the HST Cluster SN Survey which are classified 
as early-type in §5 (we also include the host of SN SCP06E12 
for analysis in this section, as it is plausibly early-type, al- 
though we leave it as unclassified in §5). 

A rough estimate of the intrinsic rate of SNe la in an indi- 
vidual galaxy can be made using the popular A+B SN la rate 
parameterization in which the rate is the sum of a term pro- 
portional to the ongoing star-formation rate and a term pro- 
portional to the stellar mass: 



R SNia = AM + BM 



(2) 



Several choices of A and B are available in the literature. 
We choose the values derived from photometric estimates 
of SN host masses and star formation rates in the SNLS: 
A = 5.3 x lrr^yr^rvlo" 1 , B — 3.9 x 10- 4 M Q _1 (Sul- 
livan et al. 2006). The values for A and B depend on which 
initial mass function (IMF) is assumed. The values above are 
derived assuming a Kroupa (2001) IMF, which is very similar 
to the Chabrier IMF with which the BC03 SSPs we have been 
using are generated. 

The cosmic SN CC rate has been shown to be propor- 
tional to the cosmic star formation rate. By dividing the cos- 
mic SN CC rate found in Bazin et al. (2009) at z = 0.3 
(RsNcc = 1-63 x 10 _4 (/T-7o 1 Mpc) _3 ) by the cosmic star for- 
mation rate at z — 0.3 found in Hopkins & Beacom (2006) 
of O.O3M0yr -1 Mpc~ 3 (and converting from the modified 
Salpeter IMF used there to a Kroupa IMF) we can obtain a 
similar formula for SNe CC: 



RsN c 



CM 



(3) 



with C = 6.52 x 10- 3 M o _1 . 

To estimate the star formation rate in these galaxies we rely 
on the the spectroscopic [O II] luminosity. As mentioned ear- 



lier, while [O II] may indicate star formation, it may also in- 
dicate LINER activity in an otherwise passive galaxy. The 
[O II] -inferred star formation rate should thus be viewed as 
an upper limit to the true star formation rate, which may be 
much less or even zero. We investigate the [O II] -inferred star 
formation rate parameterization of Kennicutt (1998) scaled to 
a Kroupa IMF: 

SFR [0 ^(Aioyr- 1 ) = 2.65 x 10- 41 L [o „](erg s" 1 ) (4) 

The [O II] -inferred star formation rate is sensitive to galac- 
tic metallicity and dust. These are in turn correlated with the 
galaxy's stellar mass. Gilbank et al. (2010) have investigated 
the correlation of the nominal [O II] -inferred star formation 
rate with other more robust star formation rate indicators in 
different mass bins and derived an empirical correction as 
a function of galaxy stellar mass that accounts for trends in 
galaxy metallicity and dust: 



SFR corr 



SFR n 



a tanh[(a; — b)/c] + d 



(5) 



where SFR corr is the corrected star formation rate, SFR nom 
is the nominal star formation rate from Equation 4, x = 
log(M*/M Q ), a = -1.424, b = 9.827, c = 0.572, and 
d = 1.700. Since the cluster early-type SN hosts are quite 
massive, this correction is significant. We carry out our anal- 
ysis both with and without this mass correction. 

No [O II] is detected in the spectroscopy of many of the host 
galaxies in the HST Cluster SN Survey. Rather than infer that 
the star formation is simply zero for these galaxies, we use 
Bayes' theorem to derive the probability distribution of the 
true [O II] luminosity given the observed value, taking into 
account the observational uncertainties, and enforcing a prior 
that the [O II] luminosity must be positive (i.e. we perform 
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TABLE 3 

Late-type and Unclassified z > 0.9 SN Hosts 



Name 


^host 


^■cluster 


[O II] EW (A) 


Confidence Interval 


Gini Coefficient 


Asymmetry 


«775 - 2850 


CMR residual 


Notes 


SCP Cluster SN Hosts 






















SN SCP06B4 


1.116 


1 12 


—48.2 


[—60.0 


—36.0] 


n 4-0 010 


1 AA 4- 01 1 


981 4- 04^ 

U.ZO-L HI U.U'iU 


_o 


709 4_ o 060 
1 ZZ ZlZ u.uuu 




SN SCP06E12 a 


1.026 b 


I .Uzo 


Q(i A 

— Z0.4 


r a i a 
y — 4f .4, 


n m 
— U.UJ 


U.ooo ± U.U1 / 


U.UOO ± U.U1 ( 


U.yoO ± U.140 


i n 
+U 


4_ n 1 z.A 
U( t zt U.104 


T 

Vu 


SN SCP06F12 


1.110 


1 i in 
I.I 1U 


on n 
— oU.U 


L — lUo.o, 


f^o 91 

— oy.zj 


n qoo 4_ n noo 
U.oyz ± U.UZZ 


n nu 1 J— n nt A 
U.Uol ± U.U14 


n aoa -i- n nc;o 
U.4Z4 ± U.UOZ 


n 


1 (I (17(1 

04y ± u.u iu 


T K 


SCP Field SN Hosts 






















SN SCP06G3 


0.962 


i .zjy 


fin i 
— 6y.l 


r a 
[— i 0.6, 


— 0Z.4J 


n /l 7q i n flic 
U.4 / o ± U.Ulo 


n (io7 4_ n nriyi 
U.Uo i zt U.UU4 


n ccn 4_ n hoa 
U.6oy ± U.Uz4 


n 


odd 4_ n i (iq 
ZUU ± U.lUo 


T 


SN SCP06N33 


1.188 


I .UZO 


— 40.0 


r A 7 ^ 
L — 4 ( .0, 


/J "3 (11 
— 4O.0J 


n A Ofi I n (115 
U.4Z0 ± U.Ulo 


n i q 1 4_ n nn7 
U. lol ± U.UU ( 


U.6UO ± U.Uoo 


n 
— U 


/109 4.(1 1 C\K 
4UZ it U.1U0 


T V, v 


SN SCP05P1 


0.926 


1.1" 


— 49.8 


[-53.0, 


-46.3] 


0.327 ± 0.006 


0.200 ± 0.010 


0.239 ± 0.035 


— 


612 ± 0.106 


L,x 


SN SCP06T1 


1.112 


0.971 


—38.1 


[—84.8, 


—0.0] 


0.333 ± 0.021 


0.119 ± 0.042 


0.167 ± 0.262 


— 


824 ± 0.280 




SN SCP06X26 


1.44 


1.101 


—4.0 


[—13.5, 


—0.0] 


0.367 ± 0.020 


0.024 ± 0.020 


0.448 ± 0.079 


— 


417 ± 0.128 


L 


GOODS SN Hosts 






















SN2002dd 


0.950 










0.345 ± 0.016 


0.114 ± 0.018 


0.187 ± 0.037 


— 


589 ± 0.106 


s,h,x,c 


SN2002fw 


1.300 










0.365 ± 0.012 


0.113 ± 0.007 


0.647 ± 0.043 


— 


215 ± 0.109 


s,h,x,c 


SN2002fx 


1.400 










0.304 ± 0.018 


0.024 ± 0.025 


0.278 ± 0.107 


— 


550 ± 0.147 


L,h 


SN2002ki 


1.141 










0.376 ± 0.009 


0.091 ± 0.006 


0.522 ± 0.039 


-0 


513 ± 0.107 


s,h,x 


SN2003eb 


0.900 










0.396 ± 0.011 


0.150 ± 0.004 


0.267 ± 0.019 


-0 


523 ±0.102 


s,x,c 


SN2003aj 


1.307 










0.419 ± 0.021 


0.143 ± 0.007 


0.589 ± 0.031 


-0 


294 ±0.105 


L,h,x,c 


SN2003ak 


1.551 










0.469 ± 0.026 


0.104 ±0.004 


0.333 ±0.026 


-0 


684 ±0.103 


L,c 


SN2003dy 


1.340 










0.405 ±0.012 


0.279 ± 0.004 


0.492 ±0.014 


-0 


417 ±0.101 


s,h,x,c 


HST04Eag 


1.019 










0.389 ± 0.008 


0.200 ±0.003 


0.370 ±0.013 


-0 


585 ±0.101 


s,h,x,c 


HST04Mcg 


1.357 










0.324 ±0.007 


0.243 ± 0.009 


0.773 ± 0.042 


-0 


124 ±0.108 


s,h,x,c 


HST04Pat 


0.970 










0.358 ±0.007 


0.140 ± 0.007 


0.584 ± 0.025 


-0 


330 ±0.103 


s,h,c 


HST04Omb 


0.975 










0.423 ±0.016 


0.287 ± 0.002 


0.191 ± 0.008 


-0 


704 ±0.100 


s,h,x,c 


HST05Fer 


1.020 










0.477 ± 0.030 


0.049 ± 0.008 


0.496 ±0.028 


-0 


410 ±0.104 


s,h,x,c 


HST05Gab 


1.120 










0.252 ± 0.019 


-0.088 ± 0.048 


0.618 ±0.187 


-0 


356 ±0.212 


s,h,x,c 


HST05Koe 


1.230 










0.284 ± 0.007 


0.119 ±0.011 


0.625 ±0.117 


-0 


318 ±0.154 


s,x,c 


HST05Red 


1.189 










0.404 ± 0.029 


0.079 ± 0.004 


0.247 ±0.014 


-0 


730 ±0.101 


L,h,x,c 


HST05Str 


1.027 










0.336 ±0.006 


0.121 ± 0.007 


0.496 ± 0.038 


-0 


457 ± 0.107 


s,h,c 


NOTE. — The [O 11 


] EW reported is the median of the posterior probability of the true EW given the observed EW, the spectroscopic uncertainties, and 


a prior that the EW be positive. 



The reported confidence interval is the smallest interval containing 90% of the posterior probability. The CMR residual is computed using the fixed-slope fits to individual clusters for cluster 
member hosts and the composite red sequence fits for field hosts. The CMR residual uncertainties reported are y/ cf nt + cr^ol wnere ^int is the measured intrinsic scatter of the parent cluster 
for cluster member hosts and a conservative 0.10 mag for field hosts. 
a Unclassified galaxy. 

k Redshift undetermined; assumed to be a cluster member at z — 1.026 

Photometric cluster redshift. 
s SN is classified as Type la from spectrum (Blakeslee et al. 2003a; Riess et al. 2004, 2007; Barbary et al. 2010) 
L SN is classified as Type la from light curve (Riess et al. 2004, 2007; Barbary et al. 2010) 
^ SN light curve passes Union2.1 cuts (Suzuki et al. 201 1) 
x SN has SALT2 XI uncertainty less than 1.0 (Suzuki et al. 201 1) 
c SN has SALT2 color uncertainty less than 0. 1 (Suzuki et al. 201 1) 



a likelihood analysis). The median star formation rates (both 
nominal and mass-corrected) from the posterior distributions 
are reported in Table 4. 

To estimate the host galaxy stellar mass, we fit BC03 SSP 
templates as in Equation 1 to the observed 1775 and z§5o mag- 
nitudes, assuming Zf orm = 3.0. The best fitting value of 
M ga i is then the initial stellar mass of the template, but the 
current stellar mass, factoring in the rapid deaths of massive 
stars, is also supplied by BC03. With only two photometric 
bands available, we only attempt this mass measurement for 
galaxies which are classified as early-types in §5. The more 
complicated SEDs of late-type galaxies will generally require 
information from more photometric bands and a more sophis- 
ticated approach to mass measurements. The statistical uncer- 
tainty of this mass measurement is limited by the uncertainty 
in our photometric measurements, both in zsso magnitude and 
in *775 — z 850 color. To this uncertainty, we also add a sys- 
tematic uncertainty (in quadrature) to capture our ignorance of 
the precise star-formation histories and metallicities of these 
galaxies. By fitting BC03 composite stellar spectra of a va- 
riety of exponential and delayed exponential star-formation 
histories and various metallicities to the 1775 and zgso pho- 
tometry, we estimate that this systematic uncertainty in mass 
is around 25% or 0.1 dex. The host galaxy masses are re- 



ported in Table 4. 

With the masses and star formation rates derived above, we 
apply Eqns. 2 and 3 to estimate the intrinsic rates of SNe la 
and SNe CC in each host galaxy. One interesting consequence 
of Eqn. 2 is that the expected increase in the SN la rate due 
to recent star-formation (the B component) is only ~ 1% - 
5% for all of our early-type hosts except that of SN SCP06U4 
for which the increase is ~ 30%. The true enhancement may, 
in fact, be smaller if the observed [O II] luminosity contains 
LINER emission. This implies that our early-type hosted 
SNe la probably have old progenitors and hence are not as- 
sociated with dust from recent star formation. 

The intrinsic rates are a start, however, since we know a pri- 
ori that the particular galaxies we are analyzing each hosted a 
detectable SN, it is more informative to estimate the apparent 
SN rates. For example, the apparent rate of a SN of type X 
given its intrinsic rate can be written: 



□apparent ^intrinsic /■detectable 

tt SN x — U SN X JSNx 



l SN x 



(6) 



where /|^° ctablc is the fraction of detectable SNe of type X 
for a particular galaxy. The zgso detection threshold for a 
point source in a single epoch of the HST Cluster SN Survey 
is about 25.3 (Barbary et al. 2010). However, since the SNe in 
question have passed quality cuts requiring detection in mul- 
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tiple epochs and multiple bands we adopt an effective zg50 de- 
tection threshold for SN peak brightness of ~ 24.7. For each 
host galaxy, we estimate the fraction of SNe la and SNe CC 
with peak zgso magnitude brighter than this detection thresh- 
old using the rest frame R -band luminosity functions derived 
from the Lick Observatory Supernova Search (LOSS) Li et al. 
(201 1). Using the maximum light spectral templates of Hsiao 
et al. (2007) for normal SNe la and Nugent et al. (2002) for 
SN1991bg-like and SN1991T-like SNe la, SNe Ibc, SNe II- 
N, SNe II-P, and SNe II-L we ^-correct our 2:350 detection 
threshold to the rest frame /?-band for each SN subtype. We 
estimate the fractions of detectable SNe la and SNe CC as 
the fractions of LOSS SNe la and SNe CC that exceed these 
^-corrected thresholds: 



detectable 



JSN_ 



Ni Nx /N SNx 



(7) 



where Ng N and N$n x Slle me number of SNe of type X 
brighter than their respective i?-band threshold and the total 
number of SNe of type X, respectively. We find that even 
at z = 0.9, none of the 25 SNe Ibc and only four (three 
SNe II-N and one SN II-L) out of 80 SNe II from LOSS would 
have been detected in our survey. On the other hand, our sur- 
vey would have discovered 55 out of 74 LOSS SNe la, or 
53 out of 59 if sub-luminous SN1991bg-like and SN2002cx- 
like SNe la are excluded. We again employ Bayes' theorem 
to derive the probability distribution of the expected value of 



N. 



SN X 



(and hence the probability distribution of / 



detectable 
SN X 



) 



assuming the observed number is Poisson distributed and a 
positive flat prior. 

With probability distribution functions for the masses (as- 
sumed to be log-normally distributed), [O II] luminosities, 
and detectable fractions of SNe la and SNe CC for each 
galaxy, we can use a Monte Carlo simulation to sample from 
these distributions and derive a probability distribution func- 
tion of the apparent rates of SNe la and SNe CC. Given the 
apparent rates, the probability that a SN is a SN CC is: 



P(CC) 



^apparent 
n SN cc 



^apparent 
U SN CC 



^apparent 
U SN Ia 



(8) 



The results of this analysis are shown in Table 4. Even un- 
der the assumption that [O II] traces star formation and not 
LINER activity, five of eleven SNe have P(CC) < 0.02, 
and four others have P(CC) < 0.15. One of the two re- 
maining SNe, SN SCP06U4, is a spectroscopically confirmed 
Type la. The final remaining SN, SN SCP06E12, was already 
considered uncertain due to its lack of spectroscopic redshift, 
though its light curve is consistent with a Type la at the red- 
shift of the cluster in the same field of view. In fact, for almost 
all of the early-type hosted SCP SNe, the already strong type 
constraints derived here from just the host galaxy information 
can be supplemented with additional constraints by consider- 
ing the full SN light curve shapes, colors and magnitudes. Ten 
out of the eleven early-type hosted SCP SNe analyzed in this 
section have evidence independent of their hosts or spectra in- 
dicating that they are SNe la (for details, see paper II of this 
series: Barbary et al. 2010). (These SNe are indicated with a 
"L" in the Notes column of Table 4). The one remaining SN, 
SN SCP06K18, does not have sufficient early time coverage 
to constrain its type through its light curve. However, as indi- 
cated in Table 4, its massive early-type host and lack of [O II] 
emission strongly suggest it is a Type la. 



7. EARLY-TYPE HOST DUST CONSTRAINTS 

In this section we use measurements of the red sequence 
scatter to place constraints on the scatter of reddening affect- 
ing these galaxies. Under the assumption that at least some 
red-sequence galaxies have very low dust content, this scatter 
constraint also sets an absolute scale for dust. 



7.1. Red-sequence scatter 

Color-magnitude diagrams of each of our clusters are pre- 
sented in Figure 15 in Appendix A. For most clusters, the 
red sequence is clearly visible as an overdensity of morpho- 
logically early-type galaxies with 1775 — z§5o ~ 1- For 23 
of the 25 clusters, the marginal likelihood distributions for 
the cluster CMR intrinsic scatters show clear maxima be- 
tween 0.0 and 0.1 magnitudes. The intrinsic scatter marginal 
likelihood distributions for clusters ISCS J1434.7+3519 and 
ISCS J1433.8+3325 are essentially flat from 0.0 up to 0.3 
magnitudes (the largest scatter for which we sampled the pos- 
terior likelihood) because the clusters contain too few red 
early-type galaxies to constrain the fit. Using the CMR fits 
as a baseline with which to combine clusters, we find that the 
intrinsic scatter of the stacked "color-magnitude" diagram is 
0.046 mag (Figure 9). Turning to individual clusters, we find 
that the best-fit intrinsic scatters of three clusters are consis- 
tent with zero. The best-fit intrinsic scatters of the remaining 
clusters range from 0.023 mag to 0.087 mag (see Table 1 for 
the best-fit values and uncertainties of each cluster's intrinsic 
scatter). 

To compare these observer-frame CMR scatter measure- 
ments to published scatter measurements at other redshifts, 
we first convert to the rest-frame a(U — V) z= o scatter. Us- 
ing BC03 SSPs, we construct a library of mock galaxy spectra 
with a variety of formation redshifts, exponentially declining 
star-formation timescales, and metallicities. For each cluster 
redshift, we determine the relation between observer-frame 
*775 — z &50 color and rest-frame U — V color through syn- 
thetic photometry of these library spectra. The slope of this 
relation evaluated at the observed CMR color ranges from 
~ 1.9 to ~ 3.2 depending on the redshift; this is the de- 
sired multiplicative conversion factor from (7(1775 — Z85o) to 
a(U — V) Z =Q- Our rest-frame a(U — V) z —o estimates lie in 
the range ~ 0.02 — 0.19 mag with a typical value of ~ 0.10 
mag and are listed in Table 1 . This range and typical value are 
consistent with the compilation of rest-frame cluster scatters 
presented in Jaffe et al. (2011), which span the redshift range 
0.0 < z < 1.46. In the rest of this section we use the proper- 
ties of cluster red sequences to constrain the amounts of dust 
present in cluster galaxies. 

A precision upper limit on dust along the line of sight to 
each SN la seen in the early-type hosts analyzed here is not 
possible with this particular data set (this rest-frame color, in 
particular, is not as sensitive to dust as a broader color baseline 
would be); nonetheless, it is useful to outline a rough analysis 
to see that dust limits on these z > 1 cluster-hosted SNe are 
consistent with the low dust values expected for red and dead 
hosts as outlined in §2.2. 

We can estimate such limits using the expectation that dust 
will both redden the intercept of the CMR and broaden the 
scatter of CMR residuals. While it is conceivable that one 
could construct dust distributions such that blue star-forming 
galaxies are reddened to precisely fall on the CMR (and sub- 
sequently decrease the CMR scatter), we take this to be an un- 
likely scenario requiring fine tuning. In fact, since the effects 
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TABLE 4 

A Priori Rate Estimates for Early-type SCP SN Host Galaxies 



Name 


log(Mass) 


median SFR 

Meyr" 1 


median Rg Npc '' 
10" 3 yr 


median RsNj a a 
10" 3 yr 


P{CC) 


Notes 


Nominal Star Formation Rates 
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SN SCP06K18 


11.8 


0.17 


0.01 
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<U.U1 
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10.4 
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SN SCP06U4 


11.1 
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Mass-Corrected Star Formation Rates 














SN SCP05D0 


10.9 


1.93 


0.47 


( 12.44) 


3.11 


( 4.51) 


0.13 


s.L 


SN SCP05D6 


11.4 


9.26 


0.84 


( 59.57) 


9.20 


(17.35) 


0.08 


L 


SN SCP06A4 


10.7 


1.04 


0.13 


( 6.70) 


1.79 


( 2.87) 


0.06 


L 


SN SCP06C0 


11.2 


7.99 


1.51 


( 51.37) 


8.13 


(12.13) 


0.15 


L 


SN SCP06E12 b 


9.9 


1.80 


0.41 


( 11.55) 


0.77 


( 1-18) 


0.35 


L 


SN SCP06G4 


11.2 


0.42 


0.03 


( 2.71) 


4.78 


( 9.27) 


<0.01 


s.L 


SN SCP06H5 


11.6 


1.16 


0.15 


( 7-46) 


12.42 


(21.32) 


0.01 


L 


SN SCP06K0 


11.4 


1.18 


0.09 


( 7.57) 


5.09 


(14.86) 


0.02 


L 


SN SCP06K18 


11.8 


0.61 


0.05 


( 3.92) 


11.14 


(32.58) 


<0.01 




SN SCP06R12 


10.4 


0.07 


<0.01 


( 0.44) 


0.81 


( 1.38) 


0.01 


L 


SN SCP06U4 


11.1 


27.01 


6.87 


(173.73) 


11.94 


(17.90) 


0.36 


s.L 



a Rates in parentheses are intrinsic quantities. Rates not in parentheses are apparent and factor in SN detectability. 



Redshift undetermined; assumed to be a cluster member at z — 1.026. This galaxy is left unclassified in §6.2. 
8 Spectroscopically confirmed SN la. 

^ Light curve shape, color and magnitude consistent with Type la (Barbary et al. 2010). 



of dust and age show different redshift evolution (see Figure 
10), the amount of dust required in such a scenario would 
need to be fine-tuned to the specific cluster redshift. Further- 
more, the typical timescale for dust destruction, ~ 10 7 — 10 8 
yr (Jones et al. 1996; Temi et al. 2007; Clemens et al. 2010), is 
significantly shorter than the typical ages of even the youngest 
z ~ 1.2 red-sequence ellipticals (~ 10 9 yr)(Rettura et al. 
2010), so we do not expect a difference in dust between blue 
and red extremes of the red sequence to play a role in the 
CMR residual scatter. Finally, a dust distribution in which 
each galaxy is reddened the same amount is a similarly un- 
likely fine-tuning scenario, and is inconsistent with findings 
in nearby ellipticals, in which infrared-inferred dust masses 
vary by several orders of magnitude (Tran et al. 2001; Temi 
et al. 2004, 2007). 

7.2. Dust scatter 

We begin by considering the dust constraints obtainable 
from the CMR scatter and then return below to the constraints 
from the CMR intercept. Although both age and metallicity 
affect the unextinguished colors of galaxies, following the lit- 
erature (Bower et al. 1992; van Dokkum et al. 1998; Blakeslee 
et al. 2003a; Mei et al. 2009; Jaffe et al. 2011), we assume 
that age is the dominant variable affecting the unextinguished 
CMR residual scatter (and that metallicity is primarily respon- 
sible for the CMR slope). We therefore assume a constant 
metallicity model in which differences in age and differences 
in dust reddening create the observed CMR scatter. 

As mentioned above, the observer-frame CMR scatter 
c(«775 — Z85o) attributable to differences in galaxy ages 
evolves with redshift (Figure 10) due to the shifting overlap 
of filters with spectral features. In particular, above z ~ 1.2, 
the 1775 and zsso filters no longer bracket the redshifted age- 
sensitive 4000A Balmer break and the predicted scatter for a 



given age distribution consequently sharply decreases. On the 
other hand, because the Cardelli et al. (1989) reddening law is 
smooth over the relevant wavelengths, the scatter attributable 
to a given distribution of dust is nearly constant with redshift. 
This complementary behavior permits us to simultaneously 
fit both the average cluster star-formation history and average 
cluster dust distribution using the scatter - redshift data gath- 
ered in §4.3. 

We investigate the age-scatter model first described by 
van Dokkum et al. (1998) in which early-type galaxies are 
formed in delta function starbursts at times uniformly dis- 
tributed in the interval(£(zo),£ z — idclay)- Heret(zo) indicates 
the epoch when cluster galaxies first form at redshift z , t z in- 
dicates the epoch at which the cluster is observed, and tdciay 
corrects for progenitor bias (van Dokkum & Franx 2001) by 
allowing time for galaxies which recently ceased star forma- 
tion to become red and evolve morphologically so that they 
will be identified as red-sequence early-types. At any red- 
shift, the scatter due to age variation can be computed from 
a population of synthetic BC03 galaxies generated with ages 
drawn from the distribution defined by zq and <dciay The 
amplitude of the effects of age variation may in general de- 
pend on the assumed metallicity. To test this dependency, we 
have carried out the present analysis for both Z = 0.02 and 
Z = 0.05, which are appropriate for the masses of early-type 
galaxies studied here (Trager et al. 2000). The final results are 
essentially the same with either choice of metallicity. To addi- 
tionally account for the effects of dust, we use the parameter 
a(E(B — V) z= o), the rest-frame scatter of B — V reddening. 
This is related to a(E(iT!^ — z^so)) by a redshift dependent 
multiplier, which we compute through synthetic photometry 
of typical early-type galaxy SEDs with and without dust, sim- 
ilar to the analysis above in which we compare our scatter 
measurements to the literature. The measured cluster scatters 
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m residual 



FIG. 9. — Stacked "Color-Magnitude" Diagram. The x-axis is galaxy £850 — m* (where the characteristic magnitude, m*, is modeled to passively evolve 
with redshift as described in the text) while the y-axis is galaxy color-magnitude relation residual: (1775 — 2350) — CMR. Open grey circles indicate galaxies 
which do not pass our broader morphological cuts; solid black circles do. The solid red line indicates zero CMR residual. The dashed red lines indicate the lcr 
intrinsic scatter of the stacked red sequence CMR. The red histogram on the right is the projection of CMR residuals for the stacked cluster galaxies. The blue 
histogram is the projection of CMR residuals for the GOODS control fields. The 2 > 0.9 SN la hosts brighter than m* + 1.5 are overplotted as stars (HST 
Cluster SN Survey) or bowties (GOODS field SN surveys). Red (gold) stars and bowties indicate field (cluster) SN la hosts that meet our morphological and 
spectroscopic early-type criteria as described in §5. Blue stars and bowties indicate field or cluster SN la hosts that do not meet these criteria. Both the cluster 
and field early-type SN la hosts are consistent with the red sequence and may have similar properties such as low dust content and old stellar populations. 



are fit with three parameters - Zq, idday, and <r(E{B— V) z =o) 
- using a Markov Chain Monte Carlo. We find that the value 
of Zq has little effect on the model likelihood within the range 
4.0 < zq < 10.0. The most likely value for ideiay is 2.0 Gyr, 
which corresponds nicely to youngest spectroscopically mea- 
sured red-sequence galaxy age in cluster RDCS J 1252-2927 
of 1.9 ± 0.5 Gyr (Rettura et al. 2010). The maximum likeli- 
hood value for cr(E(B - V)*= ) is 0.031 mag, with 68% and 
95% upper limits of 0.035 mag and 0.042 mag, respectively. 

If clusters with higher intrinsic CMR residual scatter at 
a given redshift (such as ISCS J1434.4+3426 which hosts 
SN SCP06H5) in fact contain dustier galaxies, then their 
CMR intercepts should also be redder. We have further in- 
vestigated this possibility by searching for a correlation be- 
tween Act (1775 — Z85o) an d A&22 using pairs of clusters with 
I Az| < 0.02. We find that the pair members with higher scat- 
ter are as likely to have bluer intercepts as redder intercepts. 
Thus the difference between the measured intrinsic scatter and 
the age-attributed intrinsic scatter should be interpreted only 



as an upper limit to the dust content of a cluster, and not as a 
dust detection. 

7.3. Absolute scale for dust 

So far we have estimated approximately how much CMR 
residual scatter is attributable to dust. However, for SN mea- 
surements we are principally interested in the corresponding 
amount of reddening or extinction. To convert from residual 
scatter to reddening, we consider a toy model for dust. In a 
mid-infrared study of the cluster RXJ 17 16.4+6708, Koyama 
et al. (2008) found that only ~ 10% of red cluster galaxies 
showed emission at 15/xm from dust (compared to ~ 50% 
of blue cluster galaxies). As discussed in §2.2, studies of 
nearby early-type galaxies also suggest that at least some have 
very small levels of extinction. We therefore assume that 
most cluster early-type galaxies are unaffected by dust and 
that the distribution of dust in the remaining galaxies follows 
a half-Gaussian (we have also considered exponential, tophat 
and Dirac delta function distributions for the remaining galax- 
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FIG. 10. — Measured early-type red-sequence scatter and fit assuming that 
CMR scatter is produced by scatter in galaxy ages and dust. The data points 
and error bars are our measurements and confidence intervals for the 23 clus- 
ters with scatter constraints. The blue line indicates the scatter attributable 
to age differences in our model, the red line indicates the scatter attributable 
to dust, and the magenta line indicates the combined scatter that is fit to the 
data. Clusters which hosted SNe la in early-type galaxies are circled. The 
steep evolution of the age-related scatter at 2 ~ 1.15 is due to the shifting 
of the 400()A Balmer break through the z$so filter. The SEDs of passive 
galaxies of different ages have smaller 4775 — zg^o color differences above 
this redshift because the filters no longer bracket the break. In contrast, the 
smooth dust extinction law produces scatter roughly independent of redshift. 

ies; the differences are small). This highly simplified model 
captures the main features we wish to study here: a fraction 
of dusty galaxies (and hence a remaining fraction of nearly 
dust-free galaxies) and a characteristic scale for the amount 
of dust. The full probability distribution function for our toy 
dust model is: 



/(X;/ d)W ) = (l 



f d )6(X) + f d 2 ^ J ^exp 



-X 2 



(9) 



where the parameters and fid specify the fraction of dusty 
galaxies and the mean dust content of the subset of dusty 
galaxies, respectively; S(-) is the Dirac delta function, which 
represents the dust-free galaxies (this term ensures that the in- 
tegral over the full distribution is one), 9() is the Heaviside 
step function, and X — E{B — V) z —o quantifies the rest- 
frame galaxy reddening due to dust. 

The scatter of our toy model, which can be calculated ana- 
lytically, is a(X) = jt/dVTl - fd)fd (Figure 11). Thus, for 
<j(X) = 0.035 (the 68% upper limit) we obtain constraints 
of h < 15% if \i d is large (> 0.08) or fi d < 0.06 if f d 
is large (> 0.3). The implication for SNe hosted by red- 
sequence galaxies is that either the probability of being af- 
fected by dust at all is small or alternatively that the amount 
of dust affecting the SN is small. Turning our attention to in- 
dividual SN host clusters, we see that the scatters of clusters 
RCS J234526-3632.6 (hosts SN SCP06U4), XLSS J0223.0- 
0436 (hosts SN SCP06R12), and ISCS J1438. 1+3414 (hosts 
SCP06K0 and SCP06K18) at redshifts 1.037, 1.215, and 
1.414 respectively, are below or just above the age-related 
scatter curve in Figure 10. Cluster RCS J022 144-032 1.7 (host 
of SCP06D0) at redshift 1.017 has more scatter than the three 



FIG. 1 1 . — Contours of constant CMR residual scatter attributable to dust, 
(j{X), in our toy model, in magnitudes. The model parameters f d and fid 
specify the fraction of dusty galaxies and the mean dust extinction of the 
subset of dusty galaxies, respectively. The solid and dashed contours indicate 
the la and 2cr upper limits on cr(X) respectively. The fid - fd parameter 
space above these curves are excluded at the indicated level, implying that 
either the amount of dust in an individual cluster member is low (fid is small), 
or that the probability that an individual cluster member is dusty is small (fd 
is small). 



other clusters within \ Az < 0.02 1 , but has the bluest intercept. 
Similarly, cluster ISCS J1434.4+3426 (host of SCP06H5) at 
redshift 1.241 has large intrinsic scatter but no more red an 
intercept than cluster RDCS J1252.9-2927 at redshift 1.237, 
which has only half as much scatter. 

The dust constraints derived to this point apply to cluster 
early-type galaxies. However, Figure 9 demonstrates that the 
field early-type SN hosts have similar red-sequence residu- 
als as cluster early-type SN hosts. The weighted mean resid- 
ual of cluster early-type SN hosts is 0.03 ± 0.03 whereas 
the weighted mean residual of field early-type SN hosts is 
—0.01 ± 0.02. Thus, the field early-types are actually slightly 
bluer than cluster early-types and unlikely to have significant 
reddening. 

8. SN IA CORRELATIONS WITH HOST GALAXIES 

In this section we look at correlations between the prop- 
erties of SN la light curves and the properties of SN la host 
galaxies at high redshift. 

8.1. Stretch and color 

One easily observed distance independent property of 
SNe la is the rate at which their light curves rise and fall. This 
property has been parameterized several ways in the litera- 
ture: as the stretch, s, which linearly stretches the time axis of 
a template light curve to match the observed (rest-frame) light 
curve (Perlmutter et al. 1997), Ato^B), which is the decline 
in S-band magnitudes from peak to 15 rest-frame days later 
(Phillips 1993), or x\, which is roughly the coefficient of the 
first component in a principle component analysis of SN la 
spectral time series (Guy et al. 2007). The light curve shape 
has been shown to correlate with SN la host type at lower 
redshifts with the light curves of SNe la hosted by early-type 
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galaxies rising and falling more quickly than the light curves 
of SNe la hosted by late-type galaxies (i.e. early-type galax- 
ies host SNe la with smaller stretch) (Hamuy et al. 1996; Gal- 
lagher et al. 2005; Sullivan et al. 2006). Here we compare 
SN stretch as a function of host type for our z > 0.9 SN la 
dataset. 

In Suzuki et al. (2011) we use SALT2 to fit light curves 
of z > 0.9 SNe la from the HST Cluster SN Survey and the 
GOODS SN Surveys in the ACS 2775, ^850 filters and also the 
Fl 10W and F160W filters of NICMOS. To convert the SALT2 
light curve shape parameter xi to stretch, we use the cubic 
relation from Guy et al. (2007). The distributions of stretch 
(and xi) for SN la subsets split by host type are shown in 
Figure 12. To prevent poorly measured SNe la from influenc- 
ing our results, we have only included SNe la with x\ uncer- 
tainty less than 1.0, which is roughly the size of the histogram 
bins. The SNe la whose hosts we classify as passively evolv- 
ing early-type galaxies show smaller stretch than the SNe la 
whose hosts we classify as late-type galaxies, consistent with 
lower redshift results. A Kolmogorov-Smirnov (K-S) test re- 
veals that the probability that the stretch values from the two 
host subsets are drawn from the same distribution is < 0.01. 

We have also compared the light curve shape distributions 
of these z > 0.9 SNe la to that of the first year of SNLS 
SNe la at 0.2 < z < 0.8 presented in Sullivan et al. (2006, 
Figure 12). We assume that the individual values of stretch in 
the histogram plotted there are uniformly distributed within 
their bins. We find that the K-S probability that the s values 
of the early-type hosted SNe la analyzed in this paper and the 
s values of the passively hosted first year SNLS SNe la are 
drawn from the same distribution is quite high at 0.91. The 
K-S probability comparing the s distribution in z > 0.9 late- 
type galaxies to that from the star-forming galaxies of first 
year SNLS SNe is also suggestive of similar distributions at 
0.63. Finding similar demographics across such a wide red- 
shift baseline instills confidence that SNe la are truly stan- 
dardizable candles over this range. 

Likewise, we can compare the distribution of the color pa- 
rameter in both host subsets. In Figure 13 we compare the 
distribution of color for SN subsets of different host types. 
The K-S test reveals that the probability that the c values from 
the two host subsets are drawn from the same distribution is 
0.75, also consistent with findings at lower redshift (Lampeitl 
et al. 2010). One possible explanation for this result is that 
the dominant source of color for these SNe la is intrinsic or 
circumstellar and that host galaxy extinction plays a smaller 
role. 

8.2. SN la Hubble residual vs. host mass 

Recent analyses of SN la hosts have indicated that after 
corrections for stretch and color have been applied, SNe la 
are brighter in galaxies with more stellar mass (Kelly et al. 
2010; Sullivan et al. 2010; Lampeitl et al. 2010). Kelly et 
al. have parameterized this effect both as a step function at 
log — 10.8 and as a linear function. In §6 we argued that 
the relatively simple SEDs of early-type galaxies (as opposed 
to the generally more complicated SEDs of late-type galaxies) 
allow us to measure stellar masses of early-type hosts using 
just the ^775 and z 850 photometry collected here. Thirteen of 
our z > 0.9 SNe have both host mass measurements and Hub- 
ble residuals (uncorrected for host mass) from Suzuki et al. 
(2011). To extend this sample, we have also used the masses 
of z > 0.9 GOODS SN hosts from Thomson & Chary (201 1), 
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FIG. 12. — Distribution of SN light curve shape parameters stretch (bottom 
axis) and x\ (top axis). The x\ axis labels are computed using the cubic 
polynomial conversion from stretch given in Guy et al. (2007). Both axes 
apply to both panels. Top: Light curve shape distribution for the HST- 
observed z > 0.9 SNe la discussed in the current paper. The red (blue) 
histogram corresponds to SNe la with hosts classified as early (late) type. 
Bottom: Light curve shape distribution for SNe drawn from the first year of 
SNLS as given in Sullivan et al. (2006). The red (blue) histogram corresponds 
to SNe la with hosts classified as passively-evolving (star-forming). 
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FIG. 13. — Distribution of z > 0.9 SN la color for SNe with hosts classi- 
fied as either early-type (red) or late-type (blue). 

which are derived from SED fits to ACS optical, IRAC near 
infrared and MIPS 24 jim photometry. The Hubble residuals 
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for these SNe are also available from Suzuki et al. (201 1), in- 
creasing the sample from 13 to 23. Five galaxies with mass es- 
timates are in common between the Thomson & Chary (201 1) 
dataset and the early-type hosts analyzed here. For four of 
these galaxies, the difference in mass is less than la, and it is 
less than 2a for the fifth. In Figure 14, we plot Hubble residu- 
als against host stellar masses, using our own mass estimates 
for the galaxies with two mass estimates. A suggestion of 
the trend of negative residuals for higher mass hosts is appar- 
ent. A linear fit to these data results in a slope measurement 
of —0.098 ± 0.085 mag dex -1 , (i.e. a l.ler detection signif- 
icance). The next generation of SN la studies (such as the 
SNe la discovered and followed with Wide Field Camera 3 
as part of the HST Multi-Cycle Treasury Programs) may have 
sufficiently better data to confirm this trend at high redshift. 
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FIG. 14. — Hubble residual as a function of host galaxy stellar mass. 
Open blue squares are late-type galaxies with masses from Thomson & Chary 
(2011). Filled red circles are early-type galaxies with masses derived from 
SED fits to 1775 an d 2 850 as described in §8. The solid line is the best-fit 
linear relationship with slope —0.098 ± 0.085 mag dex -1 . The dashed line 
and step function are the relations plotted in Kelly et al. (2010) when using 
the SALT2 light curve fitter. 



9. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the relation between SNe la 
and their host galaxies at high redshift (z > 0.9). In total, 
we study 16 SNe la (seven in clusters, eight in the field, and 
one with uncertain cluster membership) from the ACS cluster 
fields and 24 SNe la from the GOODS fields (Dawson et al. 
2009; Gilliland et al. 1999; Blakeslee et al. 2003b; Riess et al. 
2004, 2007). Using color-magnitude diagrams constructed 
from 1775 and z§50 photometry and measurements of quan- 
titative morphology, we have developed a technique to seg- 
regate high-redshift, early-type galaxies from high-redshift, 
late-type galaxies in these fields. Additional measurements 
of the [O II] emission line strength in the spectra of SN host 
galaxies were used to estimate galaxy star formation rates and 
the rate of core collapse SNe. This classification scheme was 
applied to the hosts of the 40 SNe, first in conjunction with the 
[O II] emission line strengths to determine the likelihood of 
early-type hosted SNe being Type la, then to identify trends in 



the supernova properties with host galaxy type. Of the SNe la 
studied here, 17 have host galaxies classified as early-type, 22 
have host galaxies classified as late-type, and one galaxy has 
a host which we have left unclassified. 
From our analysis of these data, our principle results are: 

(i) SN Classification A review of the literature (Hakobyan 
et al. 2008) shows that the confidence level is quite high for 
the classification of a SN la based on the identification of its 
host galaxy as early-type. We here confirm, on the basis of 
host morphology and spectroscopic limits on star formation, 
that SNe hosted by early-type galaxies in our current sample 
are likely Type la. Using the A+B SN la rate parameterization 
from Sullivan et al. (2006) (Eqn. 2), and by deriving a sim- 
ilar formula for SNe CC (Eqn. 3), we estimate the intrinsic 
rate of SNe la and SNe CC in individual early-type galaxies. 
By multiplying the intrinsic rates by the detectable fractions 
(derived from the detection threshold of the survey and the 
luminosity functions of different types of SNe) we calculate 
the rate of detectable SNe la and SNe CC for each early-type 
SN host. For most of the early-type hosted SNe in our sur- 
vey, the probability of misclassification is very small, even if 
relying on only the host galaxy information. The three ex- 
amples where the purely host galaxy classification is weakest 
(P(Ia) < 90%) can be traced back to observations of [O II] 
emission. If some or all of early-type galaxy [O II] emission 
is due to LINER activity instead of star formation, then these 
misclassification probabilities may be significantly overesti- 
mated. Fortunately, in this particular case two of these three 
SNe have been spectroscopically confirmed as a Type la. Of 
course, our combined classification for cosmological analy- 
sis includes the full light curve data and thus the classsifica- 
tions become much stronger (Barbary et al. 2010; Suzuki et al. 
201 1). In the HST Cluster SN Survey, ten z > 0.9 SNe in all 
can be typed by their early-type hosts, six more than a com- 
parable blank field search would have found since this many 
were found in the clusters themselves. 

(ii) Contamination from Dust The narrow color scatter 
about cluster red sequences directly constrains the amount of 
dust extinction affecting SNe in red-sequence hosts. We mea- 
sure a typical intrinsic scatter in the observed 4775 — z$5q color 
of morphologically selected early-type galaxies of ~ 0.045 
magnitudes. While other authors have modeled this scatter 
as a distribution of galaxy ages (e.g. Jaffe et al. 2011), we 
have investigated extinction due to dust as an additional pa- 
rameter. By simultaneously fitting the average cluster galaxy 
formation history and dust content, we conclude that the rest- 
frame B — V scatter of red sequence residuals attributable to 
dust is likely < 0.04 magnitudes. By modeling this distribu- 
tion as the sum of a dust-free galaxy population (with fraction 
1 — fd) and a dusty galaxy population (with fraction fd and 
mean rest-frame E(B — V) reddening fid), we show that ei- 
ther the probability that any particular red-sequence member 
is reddened by dust is small (fd < 15%) or alternatively that 
the amount of reddening is small (p,d < 0.06). These lim- 
its are comparable to early-type galaxy reddening limits set 
at low redshift obtained by exploiting the correlation of the 
Mg2 absorption line and intrinsic B — V colors of early-type 
galaxies (Ferguson 1993; Schlegel et al. 1998). We also find 
that the clusters with higher scatter have the same CMR in- 
tercepts as clusters with low scatter, contrary to expectation if 
dust were the primary factor in producing the color scatter. 

(iii) SN Light Curve Properties The light curves of z > 
0.9 SNe la rise and fade more quickly (have smaller stretch) 
in early-type hosts than in late-type hosts. The distribution of 
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SN colors is similar in both host subsets, suggesting that SN la 
color is largely intrinsic or due to circumstellar dust and not 
dominated by host galaxy extinction. These trends are well 
known at lower redshifts (Hamuy et al. 1996; Gallagher et al. 
2005; Sullivan et al. 2006; Lampeitl et al. 2010), but have not 
been demonstrated before at z > 0.9. 

(iv) Galaxy Mass and Hubble Residuals We have ob- 
tained mass estimates of the early-type classified SN host sub- 
set by fitting passively-evolving SEDs to their 2775 and zsso 
photometry. Together with infrared-derived mass estimates of 
a subset of the late-type classified SN hosts from Thomson 
& Chary (2011), we have found a suggestion of a correla- 
tion between SN la host mass and stretch-and-color-corrected 
Hubble residual (1.1a). This effect has been seen with more 
significance in three datasets at lower redshifts (Kelly et al. 
2010; Sullivan et al. 2010; Lampeitl et al. 2010) and may be 
related to the metallicity or age of SN progenitors or to ex- 
tinction by dust, all of which correlate with galaxy mass. 

While the analysis in this paper focuses on the astrophys- 
ical properties of SNe la and their host galaxies, the results 
clearly have implications for the cosmological interpretation 
of SN observations. The evolution of host demographics with 
survey strategy or with redshift can lead to biases in cosmo- 
logical parameter fits if not accounted for. In particular, many 
low redshift surveys explicitly target massive galaxies to im- 
prove their discovery rates. The cluster galaxies targeted by 
the HST Cluster SN Survey are also more massive on average 
than those in an untargeted survey and may produce SNe la 
which lead to biases in cosmological interpretation. A de- 
tailed analysis of the cosmological impact of host dependent 
properties is discussed in Suzuki et al. (2011). In that work, 
we present several other trends in SN la parameters as a func- 
tion of host type, including the values of color-magnitude re- 
lation coefficient (3 — AMb / Ac in star- forming galaxies and 
in passive galaxies, and the residual scatter about the Hub- 
ble diagram for SNe la hosted by star-forming galaxies and 
by early-type galaxies. These studies are compared to similar 
studies at lower redshift (Sullivan et al. 2003; Jha et al. 2007; 
Sullivan et al. 2010; Lampeitl et al. 2010). 

The ability to identify early-type galaxies, and their over- 
density in clusters, provides many advantages in SN cosmol- 
ogy studies at z > 0.9. SNe hosted by early-type galaxies 



are essentially all Type la, reducing the need for potentially 
expensive confirmation spectroscopy that must be obtained 
during the narrow time window that the SN is active. Clus- 
ter early-type galaxies are minimally contaminated with dust 
(either less than 15% are dusty or the reddening is less than 
0.06 mag), and can be directly checked for dust via the cluster 
red-sequence scatter. The cluster galaxies provide additional 
potential SN la progenitors and hence yield a higher rate of 
SNe la than a comparable blank field search, especially of 
early-type hosted SNe la. Finally, the relatively simple SEDs 
expected for early-type galaxies permit SN host mass mea- 
surements with fewer photometric bands than for late-type 
hosts. Targeting clusters is therefore not only the most effi- 
cient method to discover and type SNe la at z > 0.9, but is 
likely the most efficient way to collect the host galaxy infor- 
mation necessary for the next generation of z > 0.9 SN la 
studies. 
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APPENDIX 



Here we present the color-magnitude diagrams for the 25 clusters in the HST Cluster SN Survey. 



25 




FIG. 15. — Cluster color-magnitude diagrams. Small circles indicate galaxies located within 0.65Mpc of the cluster center. Dark solid circles pass our 
broader quantitative morphology cuts, light open circles do not. The slightly larger squares and triangles indicate galaxies spectroscopically confirmed to be 
cluster members, including those that lie outside of the 0.65Mpc aperture. Red squares indicate spectroscopically confirmed members lacking significant [O II] 
emission (EW > —5A.) Blue upward pointing triangles indicate spectroscopically confirmed members with significant [O II] emission. Purple downward 
pointing triangles indicate spectroscopically confirmed members found in the literature but with unpublished [O II] EWs. Stars indicate SN hosts. Blue stars are 
foreground SN hosts, gold stars are cluster members SN hosts, and red stars are background SN hosts. The solid red line is the best-fit linear color-magnitude 
relation, and the dashed lines indicate the measured scatter (including contributions from both measurement uncertainties and the intrinsic scatter). The red 
histogram indicates CMR residuals of cluster galaxies passing our broader morphology cuts falling in the magnitude range M* — 2.2 < zgso < M* + 1.0, 
which is indicated by the long tick marks on the x-axis. The blue histogram is the solid-angle-scaled analog for the GOODS control fields. 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 




FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 




FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 




FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 



35 




FIG. 15. — Color magnitude diagrams (continued). 




FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 




FIG. 15. — Color magnitude diagrams (continued). 
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FIG. 15. — Color magnitude diagrams (continued). 



